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SECTION  7 


VACUUM  ENERGY  STORAGE 

7. 1 THEORETICAL  PROGRAM 

Theoretical  consideration  of  high-energy  electron  beam  injection  and  switching 
for  the  DNA  Small  Toroidal  Pulseline  (STP)  have  been  actively  pursued  during  the  past  year. 
In  addition,  investigation  of  the  equilibrium  of  an  unneutralized  electron  beam  using  a 
fully  relativistic  rigid  rotor  model  as  well  as  the  development  of  a set  of  relativistic 
fluid  equations  for  further  equilibrium  and  stability  calculations  are  reported.  Fully 
relativistic  equilibria  are  required  to  study  the  entire  equilibrium  parameter  space, 
including  the  regime  where  the  perpendicular  motion  is  relativistic.  In  this  domain 
the  theoretical  bounds  on  stored  charge  and  u/y  consistent  with  equilibrium  at  a given 
value  of  magnetic  guide  field  can  be  properly  explored.  To  extend  our  ability  to  cal- 
culate equilibria  and  stability  in  a realistic  toroidal  geometry  requires  the  development 
of  a set  of  fluid  equations  for  non-neutral,  relativistic  plasma.  A two-mass  realization 
of  such  a set  of  equations  is  presented  in  this  report. 

In  addition  to  the  self-consistent  equilibrium  calculations,  a computer  code  has 
been  developed  to  facilitate  the  design  of  an  energetic  electron  injector.  This  code 
calculates  the  trajectory  of  a test  electron  from  the  fully-relativistlc  equations  of  mo- 
tion in  toroidal  geometry.  The  self-consistent  equilibrium  fields  of  a core  of  back- 
ground electrons  are  included,  as  well  as  a toroidal  magnetic  guide  field,  a vertical 
magnetic  field  and  a diverter  field.  This  code  permits  a study  of  the  electron  dynamics 
in  the  torus.  The  importance  of  vertical  field  for  canceling  the  electron  drifts  due  to 
toroidal  curvature  and  grad  B will  be  reported. 

Finally,  a means  of  switching  the  torus  through  the  use  of  vertical  magnetic 
fields  will  be  presented.  A design  for  achieving  the  required  switching  time  will  be 
shown  in  detail. 
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7. 1. 1 High-Energy  Electron  Injection  for  STP 


A relativistic  electrostatic  torus  can  store  times  more  charge  or  y^  times 
more  field  energy  than  is  possible  is  a non-relativistic  system  of  the  same  dimensions 
for  a given  magnetic  field  strength.  To  realize  these  advantages,  however,  the  elec- 
trons must  be  injected  with  relativistic  energy.  The  fundamental  limit  on  stored 
charge  derives  from  the  condition  for  crossed-field  electron  drift  to  be  energetically 
possible.  In  the  presence  of  the  space-charge  electrical  field,  E^,  the  self-magnetic 

field,  B , and  the  toroidal  guide  field,  B , the  electrons  execute  a poloidal  drift  with 
9 <0 

velocity  v given  by 
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2 -1  /2 

where  B = 9 E and  y = (1  - 9 ) . The  total  kinetic  energy  of  the  electron, 

ft  o r o o 

including  the  poloidal  drift,  is  given  by 


r = 
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2 

9 
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/ 2^ 
(Vg  /C)  ] 
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from  which  we  obtain  the  condition  y^v^/c  <1.  For  Tp'^g/®  ~ 1»  kinetic  energy 
of  the  electron  is  infinite  since  this  corresponds  to  v'^  = v^v^  = c in  the  reference 
frame  where  the  axial  beam  velocity  is  zero. 

The  energy  for  the  drift  motion  is  drawn  from  the  stored  field  energy  and  pre- 
sents an  "energy  penalty"  when  y^v^/c  is  near  unity.  The  fractional  energy, 

(y  - y )/y  , Is  plotted  in  Figure  7-1  as  a function  of  y v /c.  From  the  figure  we 
' ' ' O O On 

see  that  for  v v /c  < 1/2  when  the  last  electron  is  injected  at  the  torus  wall  we  ob- 
'o  9 “ 

tain  conservative  design  criterion  for  a relativistic  injection  system.  As  the  electrons 
are  magnetically  compressed  away  from  the  wall,  the  value  of  v v /c  decreases  lin- 

O H 

early  with  radius.  The  theoretical  maximum  for  y v /c  can  be  determined  by 

O ft 
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i‘ 


studying  fully-relativistic  beam  equilibria.  Such  studies  have  been  made  and  will  be 
discussed  in  a separate  section  of  this  report. 

if  ^^note  the  maximum  allowed  values  of  y by  r*,  we  obtain  an  ex- 
pression for  the  maximum  stored  charge  in  the  device.  We  write 


c y B 

O (3 

and  use  = -2Nj^e/b  wdth  Q = 2^R  N^^e  and  r = b for  the  minor  radius  of  the  torus. 
The  resulting  expression  for  the  maximum  stored  charge  is 


Q = tj'^R  b y B 
o o 


= 0.424  rty  B (kG)  millicoulombs. 

' O (5 

The  second  line  here  uses  R = 50  cm  and  b = 8, 1 cm,  which  are  typical  parameters 

for  the  STP.  This  expression  is  plotted  in  Figure  7-2  for  = 0. 5 and  B = B 

o o 

sin  (ni/2 T).  On  this  figure  we  observe  that  for  B^  = 10  kG  we  can  inject  up  to  2. 1 mC 
with  = 1.  provided  we  inject  electrons  all  the  way  to  peak  field.  To  compress  the 
electrons  away  from  the  wall,  we  must  terminate  the  injection  prior  to  the  time  of  peak 
field,  with  a corresponding  reduction  in  injected  charge.  By  using  a relativistic  in- 
jector, y > 1,  we  increase  the  allowed  charge  by  the  factor  y . 

° o 

For  a given  value  of  injected  charge,  Q,  the  energy  stored  in  the  fields  of  the 
electron  beam  is  given  by 


2 

W = 
f 2irR 


(1  + f?  ) [ 1/4  + In  b/a  ] 


where  a = the  final  beam  radius.  We  assume  Q is  the  maximum  allowed,  i.e. , 

Q = R by^B^  sin  (irt^2T),  where  t^  is  the  time  when  the  final  electron  is  injected. 
The  final  beam  radius,  a,  is  then  given  by 


10 


or 


2 2 

a B = b sin  (•nt^2T) 


b 

a 


sin(irt^2T)  J 


1/2 


Combining  these  results,  we  obtain  an  expression  for  the  final  stored  field  energ\'  in 
terms  of  the  injection  termination  time  t^, 


^ 2 2 2 2 2 

w = - — (1+^  ) trRb  y B sin  (ffty2T)[l-2Ansin('rrt /2t)  ] 

{ 8 o o o I i 

or 

= 257. 65T)^y^B^  (kG)sin^('rt^2T)  [l-2£nsin  (Trt^2T)]  Joules 

2 

for  R = 50  cm,  b = 8. 1 cm  and  fi  ~ 1.  This  expression  is  plotted  in  Figure  7-3  for 

o ~ 

r=  0. 5.  We  see  that  the  stored  energy  increases  with  t^  as  we  would  expect. 

We  now  turn  to  a calculation  of  the  power  density  available  from  switching  such 
a device.  A mechanism  for  switching  the  beam  in  one  transit  time  is  described  in  de- 
tail in  a later  section  of  this  report.  If  the  final  field  energy,  W^,  calculated  above 
is  released  in  one  transit  time,  2'ttR/c,  the  power  delivered  per  unit  area  is  given  by 


2 

(2irR/c)  (ira  ) 

Using  earlier  results,  this  expression  becomes, 

2 2 2 2 

p = — (1+|3  )y  B sln(ffty2T)  [l-2jlnsin  (ffty2T)  ] 

IGTT  000  X X 

= 1.2  X 10®T)^y^B^  (kG)sin(irt^2T)  [l-2insin  (fft^2T)  ] Watts/cm^. 
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This  expression  is  also  plotted  on  Figure  7-3  for  77  = 0.5,  and  the  result  shows  a peak 

in  output  power  density  when  injection  is  terminated  at  t^T  ~ 0, 4.  The  peak  occurs  I 

because  earlier  termination  of  injection  leads  to  a more  tightly-pinched  beam,  but 

also  results  in  less  stored  energy.  We  should  emphasize  that  these  results  are  valid 

as  long  as  the  drift  limit,  v v /c  < t),  is  the  most  stringent  limit  on  charge  injection. 

o P ■" 

We  now  examine  Figures  7-2  and  7-3  to  determine  a set  of  design  parameters 
for  a relativistic  STP.  Assume  that  we  terminate  injection  at  the  maximum  output  j 

power  point,  t^T  ~ 0.4,  and  use  = 10  kG  (the  STP  field  coils  can  reach  12-15  kG). 

The  non-relativistic  thermionic  injector  is  then  limited  to  1.  2 mC  stored  charge  cor- 
responding to  4.5  kJ  stored  field  energy'.  With  a relativistic  injector  consisting  of  a ' 

low-current  beam  of  1 MeV  electrons  (y  =3),  the  stored  charge  will  be  limited  to 

o 

3.6  mC  and  40  kJ  will  be  stored  in  the  self-fields  of  the  beam.  The  final  beam  radius 
will  be  6. 2 cm. 

The  remainder  of  this  section  will  present  some  simple  estimates  of  the  ty'pe 
of  injection  system  required.  Detailed  studies  of  the  relativistic  injector  are  under  i 

way  to  examine  the  electron  trajectories  for  various  injector  geometries  using  a i 

single -particle  model  for  the  injected  beam,  but  including  the  self-fields  of  electrons 
already  in  the  torus.  For  a low-current  injector,  such  a model  yields  a reliable  cal- 
culation of  the  injected  electron  trajectories.  The  estimates  to  be  presented  below , 
therefore,  should  be  regarded  as  an  order-of-magnitude  survey  of  the  injector  require- 
ments. 

Consider  an  injected  beam  of  1 MeV  electrons  with  current  i and  pulse  duration 

T . The  beam  radius  is  assumed  to  be  approximately  1 cm.  The  inductance  of  such 
P 

a beam  of  length  217  R is  of  order  L ~ 1 pH.  Since  the  beam  transit  time  around  the 
torus  is  2itR/c  ~ 10  ns,  the  rate  at  which  the  total  current  I within  the  torus  increases 
is  approximately 

dt  — 2irR/c 


A 
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Clearly,  a multi-turn  injection  is  mandatory.  The  induced  emf,  L(dl/dt),  associated 
with  injection  of  the  Ijeam  mast  be  balanced  by  the  transformer.  Without  a transformer, 
the  inciuced  emf  would  stop  the  beam  before  injection  could  be  completed.  If  the  injec- 
tion current,  i,  is  to.,  iante,  the  induced  emf  will  not  be  completely  cancelled  by  the 
transformer,  an.l  the  la'an.  would  l)e  degraded  in  energ\’  during  injection.  The  STP 
transformer  can  priMlui  e - 3 kV  around  the  torus.  We  therefore  have 


1.1^  ^ LiAnill^ 

■ lO-^ec 


~ 3 X 10  Volts 


i ~ 30  Amps 


for  the  average  injection  current.  From  our  earlier  considerations,  we  wish  to  in- 
ject Q ~ 3.6  mC  into  the  torus.  Since  this  charge  is  provided  by  the  injector,  we  must 


require 


Q = i T 


which  determines  the  injection  pulse  duration, 


T ~ 120  us. 

P “ 

The  magnetic  field  typically  rises  in  t ~ 3.5  ms,  and  for  injection  to  termin- 
ate at  t^T~  0.4,  we  have  t^  ~ 1,4  ms,  which  allows  ample  room  for  an  injection  win- 
dow. A lower  injection  current  with  a longer  pulse  duration  may  be  preferable  for 
obtaining  uniformity  in  the  final  beam,  and  would  be  compatible  with  the  transformer. 

7. 1. 2 Single -Particle  Orbit  Calculations 

A computer  code  has  been  developed  to  facilitate  the  design  of  a relativistic 
injection  system  for  an  electrostatic  torus.  The  trajectory  of  a test  electron  is  com- 
puted using  the  fully  relativistic  equations  of  motion  in  toroidal  geometry.  A thermionic 
(or  other)  Injector  is  assumed  to  have  supplied  a core  of  charge  within  the  torus  at  the 
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tlnu-  thf  relativistic  electron  Is  Injected.  This  core  provides  a rotational  transform 


in  the  particle  orbits  to  cancel  the  vertical  drift  due  to  the  centrifugal  and  grad  B 


forces.  Vertical  magnetic  field  can  also  be  utilized  to  cancel  the  vertical  drifts, 


thereby  eliminating  the  need  for  rotational  transform.  The  fields  due  to  the  core 


electn)ns  are  computed  from  self-consistent  equilibrium  solutions.  The  electrostatic 


and  self-magnetic  fields  of  the  core  electrons,  including  the  diamagnetic  contribution 
to  the  toroidal  field  are  all  present  in  the  equations  of  motion.  In  addition,  the  code 
provides  for  a loop  diverter  field  to  distort  the  magnetic  field  lines  in  the  vicinity  of 


the  injector,  and  a vertical  field  to  adjust  the  equilibrium  location. 


The  following  parameters  are  input  variables  in  the  code,  and  are  therefore 


arbitrary: 


1.  Test  Electron 


initial  position  of  the  test  electron 


initial  energy  of  the  test  electron 


direction  of  the  initial  momentum  vector  of  the  test  electron 


Electron  Cloud 


quantity  of  previously  Injected  charge 


radius  of  the  electron  cloud 


relativistic  factor  (y  ) of  the  electron  cloud 
' o 


3.  Geometry  of  the  Toroidal  Vessel 


major  radius  and  minor  radius 


4.  Diverter  Loop  Quantities 


current  in  loop 


radius  of  loop 


location  of  loop  center 


direction  of  loop  axis 


5.  Guide  Field,  B sin  (fft/2T) 
o 


and  t 


6.  Initial  and  Final  Time  in  Computer  Run  and  Time  Step  for  Integration 
of  Equations  of  Motion 

7.  Vertical  Field,  B 


Coordinates 


Figure  7-4  shows  the  toroidal  coordinate  system  as  well  as  local  and  global 
cartesian  coordinates.  The  problem  is  solved  in  the  global  cartesian  system  in  order 
to  facilitate  the  inclusion  of  the  diverter  field,  which  is  non-toroidal.  The  transforma- 
tion from  toroidal  coordinates  (r,  9,  o)  to  global  cartesian  coordinates  (x,  y,  z)  is 


given  by 


X = R cos  (0 
y = -R  sin  o 
z = r sin  0. 


where  R = R + r cos  0.  The  inverse  transformation  is 
o 

R = v/x  +y 


tan  o = -y/x 

/ 2 t 

r = ,/(R  - R ) + z 
* o 

tan  0 = z/(R  - R ). 

o 


The  unit  vectors  in  the  global  cartesian  system  are 


e cos  0 cos  (0  sin  0 cos  (d  - ^ sin  o 

r 0 o 


cos  0 sin  (0  + ^ sin  0 sin  o " ^ cos  o 
r 0 (0 

^ sir.  0 + ^ cos  0. 

r 0 


Given  any  vector,  A=^  A A^  + ^ A , in  toroidal  coordinates,  the  components 

— r r 9 6 o o 

in  global  cartesian  coordinates  are 


X X 


A = b'  • A = (A  cos  0 - A sin  0)  COSO"  A sin  o 


o 


r ” 0 

A = ^ • A = -A  cos  (0  - (A  cos  0 - A sin  6 ) sin  o 

y y - o r 6 

A = ^ • A = A sin0+A„cos0. 


z z 


Fields 

The  code  is  written  entirely  in  the  cgs-Gaussian  electrostatic  system  of  units. 
The  self  fields  of  the  electron  cloud  are 


E = 


r itR  a 
o 


r/a 


a/r 


r s a 


a s r ? b 


= 8 E . 
0 or 


The  toroidal  guide  field  with  the  diamagnetic  contribution  of  the  electron  cloud  is  given 
by 


0 


B sin  (nt/2T) 
o 


R 


1 - 


1 R a 
o 


(1  - rW) 


R 


R 


— B sin  (fft/2T) 
R O 


r 5 a 


a 5 r S b. 


In  addition,  there  is  an  inductive  electric  field,  E , due  to  the  time  dependence  of 

0 

B 


0 


^6 
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Diverter  Loop 


The  diverter  is  modeled  as  a single-turn  filamentary  loop.  Later  versions 
of  the  code  will  stack  such  loops  into  a finite-length  solenoidal  field.  The  diverter 
field  is  calculated  in  the  coordinate  system  shown  in  Figure  7-5.  The  vector^  de- 
fines the  axis  of  the  loop  as  well  as  the  direction  of  current  flow  in  the  loop  (by  the 

righthand  rule).  The  loop  radius,  a , and  the  current  in  the  loop,  I , are  input 

I t 

variables  in  the  code.  In  addition,  we  specify  the  location  of  the  center  of  the  loop 

(r  , 6 , o ) the  coordinates  of  an  arbitrary  point  on  the  positive  half  of  the  loop 
ji  t t 

axis  (r  ,0  , to  ).  The  point  at  which  we  wish  to  know  the  field  is  (r,  0,  o)- 

ax  ax  ax 

Using  the  transformation  equations,  we  find  the  coordinates  of  these  points  in  the 
global  cartesian  system 


(r  ,9  , (0  , - (X  , y , z ) 

ax  ax  ax)  ax  ax  ax 


(r,  9,  (O)  - (X,  y,  z) 


We  define 


i'  = (x  -x,y  -y,z  -zj 
- ' ax  I '^ax  t ax  t 


p = (x-x^,  y-y^,  z-z^) 


where  £ is  the  vector  pointing  from  the  center  of  the  loop  to  the  field  point.  The  loop 
field  is  then  given  by 
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Equations  of  Motion 


The  code  solves  the  fully  relati\dstic  equations  of  motion, 


d\ 

dt 


-E- 

/ 2,22 

m y/  1 p /me 


^ = - e j ^ 1 X.  ^ 

dt  L c 

by  a numerical  integration  technique  similar  to  the  classical  Runge-Kutta  procedure. 
Output  of  Code 

The  code  makes  two  plots,  a projection  of  the  electron  trajectory  in  the  x-  y 
plane  and  a projection  in  the  local  x-z  plane.  These  two  plots  provide  the  complete 
shape  of  the  orbit.  In  aikiition,  the  code  prints  out  the  position  and  momentum  com- 
ponents of  the  test  electron  at  each  time  step,  thereby  allowing  the  time  history  of  the 
particle  to  be  followed. 

We  have  run  a number  of  simple  problems  in  order  to  check  the  accuracy  of 

-10 

the  code.  In  one  case,  the  calculation  time  step  was  reduced  from  10  sec  to 
lo”^^  sec  wtth  no  change  in  the  output.  In  another  example,  we  reversed  the  sign  of 
the  toroidal  guide  field  and  caused  the  toroidal  curvature  and  grad  B drifts  to  change 
direction. 

One  set  of  test  calculations  were  of  pure  toroidal  drifts.  The  torus  was 

emptied,  Q = 0,  and  the  test  electron  launched  in  the  toroidal  direction  with  no  initial 

momentum  in  either  the  radial  or  the  azimuthal  directions.  The  diverter  loop  was 

not  activated,  I = 0.  In  these  calculations  we  expect  the  electron  to  execute  an  up- 
l 

ward  vertical  drift  with  speed 

« 2 3 

y /3  me 
o o . 


eR  B 

O (0 
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The  results  of  these  calculations  are  showTi  in  Figures  7-G  through  7-8.  In  each 
case,  the  distance  to  wall  is  7.746  cm,  since  the  electron  begins  2 cm  off  axis.  From 
the  known  time  of  flight  of  the  electron  to  the  wall,  we  can  calculate  the  electron  drift 
speed,  which  we  compare  with  the  theoretical  value.  The  following  table  summarizes 


the  results. 

Figure 

y 

B (kG) 

V (Theoretical) 

Z j 

(cm/ns) 

V (Computer) 
Z / 

(cm/ns) 

7-6 

2 

0 

5 

1,382 

1.383 

7-7 

3 

5 

2.456 

2.421 

7-8 

3 

10 

1.  228 

1.229 

The  small  differences  between  the  computed  and  theoretical  values  of  v^  are  entirely 
consistent  with  the  finite  Larmor  radius  effects  which  are  included  in  the  code.  The 
largest  deviation.  Figure  7-7,  corresponds  to  the  largest  Larmor  radius. 

In  the  ne.xt  set  of  test  calculations,  we  constructed  HIPAC  equilibria.  The  test 
electron  started  at  rest,  2 cm  off  the  minor  axis.  The  torus  contained  a finite  amount 
of  charge,  Q,  and  the  azimuthal  electron  drift  w'as  expected  with  speed 


V = /gr 

ft  ^ 


U! 


_S£. 


itR  a B 
o o 


Figure  7-9  shows  the  result  for  Q = 10  ;iC  and  = 1. 0001  for  the  electron  cloud. 


The  test  electron  executed  a perfect  azimuthal  drift.  The  calculated  precession  half- 
period, it/ ti!,  was  9.75  ns,  and  the  test  electron  required  9.6  ns  to  process  half  way 
around  the  minor  axis.  The  calculated  Larmor  period  and  Larmor  radius  also  agreed 
with  the  computer  results. 

Figure  7-10  shows  the  same  calculation  with  Q = 100  jiC.  The  orbits  arc  larger 
and  the  precession  faster,  as  expected.  Both  Figures  7-9  and  7-10  represent  10  ns  of 
electron  travel.  In  Figure  7-11,  we  repeat  the  calculation  with  Q = 100  but  now 
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y = 3 for  the  electron  cloud.  The  test  electron  still  begins  at  rest.  The  orbit  now 
o 

is  smaller  because  the  radial  force  is  decreased  by  the  presence  of  B . The  test  elec- 

n 

tron  drifts  in  the  toroidal  direction  due  to  x ^ B . 

9 9 

Figures  7-12  and  7-13  show  the  results  of  test  runs  in  which  the  torus  contained 
100  jiC  of  charge,  and  the  test  electron  was  launched  in  the  toroidal  direction  with  rela- 
tivistic energy.  In  Figure  7-12,  the  test  electron  was  launched  with  y ” 2,  and  the 
curvature  drift  overpowered  the  azimuthal  drift,  causing  the  test  electron  to  hit  the 
wall.  In  Figure  7-13,  we  reduced  the  test -particle  energy  to  y-  1.5,  and  obtained  a 
stable,  confined  orbit.  We  note  that  the  center  of  the  orbit  in  the  x-z  plane  does  not 

coincide  with  the  center  of  the  conducting  wall;  but  is  shifted  toward  the  inner  wall  of  , 

I 

the  torus.  The  trajectory  can  be  centered  in  the  torus  through  the  application  of  a ver- 
tical field,  B . 

z 

The  major  problem  encountered  in  relativistic  electron  injection  from  the  torus 
wall  is  that  the  electron  trajectory  lies  close  to  the  wall  and  is  not  concentric  with  the 
walls.  The  inevitable  result,  therefore,  is  that  unaided  electrons  will  hit  the  wall. 

One  exception,  which  is  not  experimentally  feasible,  is  to  launch  the  electron  from  the 
inner  wall  of  the  torus.  Since  the  trajectories  are  displaced  toward  this  wall,  an  un- 
aided electron  will  be  trapped  within  the  torus.  This  situation  is  illustrated  in  Fig- 
ure 7-14,  which  shows  the  trajectory  of  an  electron  having  y = 1.2  launched  from  the 

inner  wall  with  p = p and  p = 0.  The  electron  remains  in  the  torus  during  the 
0 o r 

entire  10  ns  duration  of  the  computer  run.  When  the  electron  is  launched  in  an  iden- 
tical manner  from  the  lower  wall  of  the  torus,  it  hits  the  upper  wall  as  shown  in 
Figure  7-15. 

To  overcome  these  difficulties,  we  use  external  magnetic  fields  to  assist  in 
electron  injection.  The  two  basic  types  of  external  field  under  consideration  are  the 
vertical  magnetic  field,  B^,  and  the  diverter  field.  The  diverter  field,  a simple  cur- 
rent loop  near  the  injection  point,  serves  to  locally  distort  the  magnetic  field  lines 
near  the  injector  so  that  electrons  can  be  introduced  along  field  lines  which  run  to  the 
interior  of  the  torus.  In  this  way,  the  electron  trajectory  is  moved  away  from  the 
wall  and  the  centering  of  the  orbit  becomes  less  critical.  The  successful  operation 
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of  the  diverter  depends  on  the  local  magnetic  perturbation  being  large  enough  to  allow 
electron  injection,  yet  small  enough  that  the  electron  is  unlikely  to  leave  the  torus 
through  the  diverter  after  only  a few  transits  of  the  torus.  The  injected  electrons  must 
remain  within  the  torus  for  ? I )is,  which  is  the  expected  duration  of  the  high-energy 
electron  pulse. 

The  vertical  magnetic  field  acts  to  displace  the  electron  orbits  toward  the 
outer  wall  of  the  torus,  thereby  allowing  us  to  make  the  orbit  concentric  with  the  torus 
w'all.  Through  the  use  of  this  field  we  can  relax  the  diverter  field,  making  successful 
containment  of  the  injected  electrons  easier. 

The  vertical  magnetic  field  B acts  by  producing  a force  on  the  electron  in  the 

z 

direction  of  the  major  radius. 


= - V B , 

R c o z 


which,  when  coupled  to  B , produces  a vertical  drift; 

(3 


B 


V = V — 

z (3  B 


(3 


This  drift  competes  with  the  vertical  drift  due  to  grad  B and  centrifugal  forces.  To 

displace  the  beam  toward  the  outer  wall  of  the  torus,  the  vertical  field  must  produce 

a V drift  which  tends  to  cancel  the  toroidal  geometry  drifts.  Neglecting  v /v 
z ‘ 1 C3 

relative  to  unity,  the  curvature  drifts  are  given  by 


zc 


3 

g ymc 
eRB 

(3 


For  these  two  driffe  to  cancel,  we  require 


B = - 
z 


ygmc 

eR 


~ - 34  •yfl  Gauss 


for  R = 50  cm.  Note  that  a negative  (downward)  vertical  field  is  required  to  center 
the  orbit  when  v^  is  positive.  This  analysis  indicates  that  the  vertical  magnetic  field 
can  be  employed  to  remove  the  need  for  a rotational  transform  by  canceling  the  verti- 
cal drifts  due  to  toroidal  geometry  (curvature  and  grad  B). 

In  Figure  7-7  we  displayed  the  upward  vertical  drift  of  a test  electron  with 
y = 3.  From  the  above  discussion,  it  appears  that  = -96  Gauss  should  just  cancel 
this  drift.  Figure  7-16  shows  that  the  expected  result  is  realized  when  we  rerun  this 
example  with  the  inclusion  of  vertical  field. 

In  Figure  7-12,  we  saw  that  Q = 100  pC  did  not  provide  sufficient  rotational 
transform  to  contain  a test  electron  with  y = 2.  By  adding  a vertical  magnetic  field, 
®z  Gauss,  however,  we  can  eliminate  the  need  for  the  rotational  transform  and 

the  test  electron  is  contained  as  shown  in  Figure  7-17.  Note  also  that  the  contained 
orbit  is  concentric  with  the  torus  wall.  The  inward  shift  due  to  the  superposition  of 
vertical  and  poloidal  drifts  has  been  eliminated  by  the  inclusion  of  the  vertical  mag- 
netic field. 

Studies  are  underway  to  inject  relativistic  electrons  from  near  the  torus  wall. 
Using  vertical  field  to  center  the  trajectory,  the  results  shown  in  Figure  7-18  w’ere 
obtained.  The  electron  was  injected  with  y = 1.5,  p = 0,  p = p from  an  internal 

r 9 o 

injector  1 cm  from  the  torus  wall.  Multi-turn  confinement  of  the  electron  is  evident 
in  the  figure,  which  show's  only  the  first  250  ns  of  the  trajectory-.  The  injector  was 
treated  as  an  obstacle  in  the  code  and  enclosed  a 1 cm  x 1 cm  x 1 cm  volume  near  the 
torus  wall.  In  this  particular  example,  the  test  electron  struck  the  injector  after 
908  ns.  Smaller  Injectors,  wiiu  electrons  emitted  close  to  the  edge  of  the  injector 
should  be  investigated  in  an  effort  o determine  if  such  an  internal  injector  can  form 
the  basis  of  a high-energy  electron  source. 

7.1.3  Self-Consistent  Relativistic  Equilibria 

There  are  several  approximate  calculations  for  self-consistent  relativistic 
equilibria.  The  simplest  approach  is  to  assume  that  the  electron  motion  in  the  "beam 
frame"  is  non-relativistlc  and  to  then  transform  the  resulting  equilibrium  to  the 


laboratory  frame  via  the  Lorentz  transformation.  The  result  so  obtained  will  be  shown 
to  be  equivalent  to  working  in  the  laboratory  frame  using  the  "two-mass"  approxima- 
tion, to  be  described  below,  in  an  expansion  to  second  order  in  the  small  parameter, 
q,  where 


This  expansion  will  be  shown  to  be  valid  when  v/y  « 1/4. 

Starting  with  the  beam  frame  (denoted  by  prime),  and  assuming  the  motion  to 
be  non-relativistic,  we  obtain  the  single-particle  constants  of  motion, 


1 2 

H'  = - mv'  - e (r* ) 

P'  = r ' m v’  - - r’  A'  (r'  ) 

6 e c e 

P'  = m v' 
z z 

where  all  quantities  are  assumed  to  depend  on  the  radial  coordinate  only.  Any  func- 
tion of  these  constants  will  satisfy  the  stationary  Vlasov  equation,  and  we  use 


where 


F(r'.  V')  = ^ °'o  ^ 


= ^“'o5  ^ mv-/-g’(r’)|6(v;) 


v’  ^ H v'  ^ + (v'  - w’r')^ 
X r 0 


g'(r’)  £ e 4>'  - A'g  + ^m  (oj'r')^  + 
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By  computing  the  moments  of  F(r  , v ),  we  see  that  this  equilibrium  corres- 
ponds to  a uniform  density  rigid-rotor  with  a finite  beam  radius,  a . The  non- 


vanishing moments  are 


n'(r’)(v’  > 


( n’  r'  < a' 

) ° ■ 

( 0 r > a’ 

I us'r'  r'  < a' 


n'(r’)<v^S  = n’(r’)  ((v'g-  u)’r’)S 


n'  g'  (r')  r'  < a' 
o 


r ’ > a' 


where  the  beam  radius  is  given  by 

g’  (a')  = 0 

together  with  the  conditions  for  self-consistency, 

g*  (d  > a' ) <0 
g (r*  < a' ) >0 

which  must  be  verified  a posteriori. 

We  may  also  define  a transverse  "thermal”  velocity 


The  self-consistent  fields  are  obtained  from  Maxwell's  equations 


1_  _d_  , d»' 

r'  dr'  dr’ 


4 IT  en’  r'  < a' 
o — 


0 r'  > a' 

4ir  en'  y'  r'/c  r'  < a 


did 

r'  A'  = 

dr'  r'  dr'  9 


0 


r'  > a 


Assuming  the  presence  of  a conducting  shell  of  radius  b',  with  (b') 
obtain 
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0,  we 


b' 


b'  . 


Returning  now  to  g'  (r’),  we  must  verify  the  conditions  for  self-consistency  and 
derive  the  relationship  between  the  beam  radius  and  cc'.  The  equation  g'(a')  = 0 implies 

2 


2 b ^ 
f ' = 2N  e i n - + 
S a 


1 .B..  (^)  - 1 (e) 

1 2 / a’  \ “ 

- me  I 

2 \ c I 


Using  this  expression,  we  obtain 


!'(r'<a')  . (‘^)  - N'e"  - J N'®"  (“t)  |‘  " ( a' ) ] 

I - ii)  T 


1 2 
- — me 
2 


which  must  satisfy  g'(r'  < a')  > 0,  or 


o_  ^ 

W'  ^ 2 


1 + 


l(“f) 


-OL 


Similarly,  for  r'  > a'  we  obtain 


g'(r'>a')  = 2N'e  in 


which  must  satisfy  g'  (r'  > a')  < 0,  or 


^ ^ 1 ^ 

to'  2 0!'  O! ' 

P P 
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Finally,  using 


^ m (v*  = g'(0)  = - N'e^  1 + 2 in  ^ + ^' 

^ X Si 


we  obtain  a quadratic  equation  for 


2 \ c 


which  has  two  solutions,  uj'^,  given  by 


2 1 + - 


c 


a'  \ <v 


2 \ c 
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The  requirement  that  be  real  introduces  an  additional  condition. 


>2  1 +-  I— ^ 
- 4 \ c 


/ o V 


2 \ a>'  a’ 

c ' P 


If  this  equation  is  satisfied  and  <v'  ) > 0,  then  a self-consistent  equilibrium  has  been 
found.  We  can  prove  this  statement  by  rewriting  the  quadratic  equation  for  u,”,  as 


. r 1 /W’  a'\  . tc’ 

“^'p  ‘^'p  *■  ^ V ® / J 2 . 


which  satisfies  the  conditions  for  equilibrium  provided  <v'  ) > 0. 

i 


2 2 
u;'  a' 

P 


and  the  condition  for  equilibrium  is 

I- 


which  we  plot  in  Figure  7-19.  From  the  fig-ure  we  see  that  q“  « 1 is  required  when 

- a/c  is  large.  Since  p/y  = .25  (w  a''c)“  is  the  ratio  of  beam  field  energy  to  par- 
'*'p  o ^ p 

tide  kinetic  energy,  we  see  that  q"  is  a small  parameter  for  equilibria  which  store 

enorgj'  mainly  in  the  beam  self  fields. 

This  equilibrium  may  also  be  computed  in  the  laboratory'  frame,  using  the 

"two-mass"  approximation,  where  the  axial  velocity  is  written  as  ~ ^ ^'z  ^o^ 

with  the  assumptions,  » v^,  | - V^|  and  assumed  to  be  constant,  in  this 

approximation  we  may  write 


d . d , 

- (ymv)  = ym  V + y m v 
dt  - “J-  "'ll 
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~ V m V + y n5  V e 
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2 

where  y = (1  - ) and  B - ^ /c. 

O O u u 


The  single -particle  equations  of  motion 


are 


y m 
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dv 
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[ dt 
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These  equations  were  derived  from  an  expansion  to  second  order  in  the  small  param- 

2 

eter,  q,  keeping  terms  to  order  q , by  Rostoker  and  Hieronymus. 

The  claim,  dA  /dt~  0,  implies  the  constancy  of  v . In  fact,  the  canonical 
z z 

axial  momentum, 


r,  3 e . 

P = y mv  - - A 
z 'o  z c z 


2 

is  the  exact  constant  of  the  motion.  To  have  dA  /dt  ~ 0 to  order  q , we  must  estab- 

z - 

lish  that  the  change  in  v during  the  particle  orbit  is  of  higher  order.  The  radial  vari- 

z 

ation  in  the  trajectory  of  a single  electron  which  executes  an  azimuthal  drift  at  velocity 
Vq  = o-'r  is 


Ar 


2v 

-1  = 


V 


To  satisfy 


^ V 


8 c 
o 


z 2 

« q 


for  second-order  constancy  in  v , we  must  require 

z 


v 1 

3 4 
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With  this  assumption,  the  single-particle  constants  of  motion  are 


1 / 2 2\ 

Hi  = - ymlv  +v„  )-e($-8A) 
2 o\r  e/  oz 


P = 


V mrv  - — r A 
'o  0 c 0 


correct  to  order  q . The  distribution  function 


y mn 

F(r.  V)  = 6 |Hi  - UJP,  - Cj6  | 


yields  the  same  equilibrium  as  was  previously  obtained  by  assuming  non- 
relativistic  motion  in  the  beam  frame. 


' 
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If  we  attempt  to  use  the  two  mass  approximation  with  the  full  constants  of 


motion. 


2 2X132 

H = — vmlv  V )+—  vmv  -eq>(r) 

2 'o  V r e / 2 'o  2 ' ' 


P = V mrv  - — rA  (r) 
e ^ f)  c e ' ' 


3 e , 

P = Y mv  - ~ A (r) 

z 'o  z c 2 


in  the  distribution  function, 


4 2 
Y m 

F(r,  V)  = -2- n 6 H - ce  P-  - C 1 6 fP  - V ^mV  1 , 

2Tr  o e^o  zo  o 


the  macroscopic  velocity  in  the  axial  direction  is  found  to  be 


<v  ) = V I 

2 O O \ Y^C 


for  r < a,  where  is  the  modified  Bessel  function  of  zero  order.  To  satisfy  the 

requirements  of  the  two-mass  approximation,  V » |v  -V  , we  must  have 

o I z ol 


or  again 


V 1 
3 <<  J • 

’'o 

With  this  inequality,  we  may  expand  the  Bessel  function  and  we  recover  the  previous 
equilibrium.  The  two-mass  approximation,  therefore,  corresponds  to  assuming  non- 

3 

relativistic  motion  in  the  beam  frame  and  appears  valid  when  i//y  « 1/4.  A more 

o 


where 


-i-  f.- 

2 1 o 

r*  ^ 


e a!'r' 


A’  + e 4>'  1 

« J 


To  carry  out  the  indicated  integrations,  the  6 -function  may  be  expressed  as 


1 2 . , , 

: ! — — 1 5(<5-o  ) 

'r'  / 2^  COSO  I o 

— y y - 1 O ! 


where 


wdth  the  requirement 


To  satisfy  this  equation,  y must  lie  in  the  range,  y_  < y £ y^.  where 


-(^y) 


L_  ( 
r-\  ^ t 


i g'(r')  ± 


as  shown  in  Figure  7-20.  The  two  intersections  y_ , y^  coalesce  to  one  when 

J 1 - (a,''r'/c)^in  which  case  y^  = - (uJ'r'/c)^  >1.  The  radius  of 

the  beam  is  defined  by 


g'(a')  = yi  - (oj'a'/c) 
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with  the  additional  requirements  that 


g'(r'  > a’)  < yjl  - (u-'r'/c) 


g'(r'  < a')  > y 1 - 


which  must  be  verified  a posteriori.  These  conditions  ensure  that  the  density  vanishes 
for  r'  > a’  and  is  finite  for  r'  < a'. 

The  moments  of  the  distribution  function  may  now'  be  evaluated  as  follows: 


n’(r’)  = K 


'm  r 

Jy 


^(a:’r7cy  - 1)  - [y  - g'  (r') 


2 ^ ■ + 


I'mc  ^ 

J 


2 [y  - g'(r’)  jdy 

I 2 2 ^ 

'(uj'r'/c)  (y  - 1)  - [y  - g'  (r')] 


which  yields  the  moments  of  F(r',  £’), 


^ 2TrK'mg'(r')  ^ g-(r-)y'^r') 

n (r  ) , 2,  3/2  g'  ^ ^ 

[1  - (a;’r'/c)  ] o 


(v’g>  = w'r' 


for  r'  < a'.  The  moments  vanish  for  r'  > a'.  We  have  defined  n'^  = n'(r'  - 0), 

/ 2,  -1/2 

g.  ^ g' (r' = 0)  and  v ' =[l-(u.’  rVc)  ] 
o i 

We  now  transform  to  the  laboratory  frame  and  write  Maxwell's  equations. 


In  the  laboratory  frame  we  have 


where 


- 1/2 


" (Vq  wr/c)  ] 


g(r)  = 


o . _ 


me 


o z 


f - t'  /y 
o ^ o 'o 


Maxwell's  equations  may  therefore  be  written 


as 


j . -4  TT  en 

7 d7  ’’  ^r  " ^ Ag  + e (4>  - «^A^)  1 ^r) 


dB^  4 TT  en 

7 

dr 


7 [C  - A + e (4>- /?  A )1  y 

^ ^ c ft  oz'^c^i 


iiUL  3 


By  multiplying  the  first  equation  by  ojr/c  and  adding  to  the  second  equation, 
obtain 


z c r 


or 


fjG 


(r)  . 


\^’e 

1 


t 


With  this  relation  we  can  reduce  the  problem  to  a single  equation  for  E , 

r 


1 E 

1 d ^ r L. 

^ ^ ^"-r  ' L " ^2 

ry  o y 

i j- 


^ 2 , 
XJ.  Y J 

\ o c 


where 


4'tTn  e 
1 _ o 

2 “ 2 

L y L. 

o ■^o  \) 

This  equation  can  be  solved  analytically  in  the  limit  y ~ 1,  where  the  perpen- 

1 ~ 

dicular  motion  is  non-relativistic.  Define  the  dimensionless  variables 

X = r/L 

o 

f = E /B. 
r 1 

tii  — te  L / C 

w E 

Po  _ P 


where 


47Tn  e 
2 o 


O.  = ^ 

1 y me 
o 


For  y ~ 1,  we  obtain 

X “ 


d f ^ 1 df 

2 X dx 
dx 


1 2 _ 

- (1  + -J  ) f = U)X 

X 


v^-ith  boundary  conditions  f (0)  = 0,  f'(0)  - -rj/2. 


The  solution  is  then 


2 - 


f(x)  = [2y^  0.'  - T)]  (X)  W' 


and  the  corresponding  potentials  are 

2 


'J>(x) 


o o 


2 - 


[1  - MX)]  + 2y^  4 


A (X)  = 8 $(x) 


= B.L^  [ I + (y/  ^ + a>  ( r - cc)  X I^IX)  ] 


2 _ 


where  <J>(0)  = 0 has  been  assumed.  From  the  equation  for  g(r)  we  sec  that 

2 


y f U!  L 

0 0 , P o V 

g = g(0)  = — ; = ( ) 
° mc“ 


and 


2 - 
2y„  03 


^ = ^o  ^ ^"4  ^ 


2 - 


y 03X  _ 2 ^ x^  _ - 2 

— [-+(toy^)  — + o.'(T7- 2 cry^  )x 


where  I (x).  I (x)  and  I (x)  are  modified  Bessel  functions, 
o 1 2 


l2(x) 


The  assumption  y ~ 1 implies  that  the  diamagnetic  contribution  to  be 
small.  In  fact,  the  solution  presented  above  is  not  exactly  self-consistent  as  may 
be  seen  by  computing  the  density  from  the  solution  for  f(x) 


n(x) 


- n 

o 


1 

TJX 


_± 

dx 


X f 


= n 

o 


[(1  - 


2 _ 

V 


I (X)  + 
o 


2 - 
U! 


V 


] 


which  yields 


I ti-i  (X)] 

o n o 

which  differs  from  the  previous  expression  for  girj/g^  by  just  the  diamagnetic  con 

tribution  to  A , and  is  consistent  with  assuming  ~ 

G 

In  terms  of 


% 


v 


or 
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the  equation  for  g(r)/gQ  becomes 


and 


Kl£l 


(1  -^)I  (X)  + 

Wo  o 


g(a) 


g 


o 


Ui 

(1  - — ) I 
w o 
o 


defines  the  radius  of  the  beam  so  that 


I (a/L  ) - g ^ 
o o q 

I (a/L  ) -1 
o o 


If  a/L  ~ w a/c  » 1,  we  obtain  w W » and  the  density  is  approximately  uniform, 
op  ~ o 

For  this  case  the  fields  are 


E (r)  ~ - 2 wen  r 
r — o 


E 

r 


B^(r) 


2Ne  w 2 

o  r 

B.  + (-) 

1 c a 


The  conditions  for  this  simple  equilibrium  are  w a/c  > 1 and  w a/c  < 1/y  , 

p o o 

which  may  be  summarized  as 

2  2 

1/4  < i//y  < (y  /r)  (fl  / w ) 

0 0 op 

where  f2  = eB  /y  me  and  B = B + 2Ne  w /c.  This  equation  implies  B > B /2, 
o o 'o  o i o i o 

which  places  a limit  on  the  diamagnetism  of  the  beam.  When  this  equation  is  satis- 
fied, the  resulting  equilibrium  is  the  same  as  that  obtained  using  the  two-mass  ap- 
proximation. This  Inequality,  therefore,  gives  the  condition  for  validity  of  the  two- 
mass  approximation. 
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We  now  turn  to  the  numerical  solution  of  the  field  equation  for  arbitrary 
In  terms  of  the  dimensipnless  variables,  we  have 


_d 

dx 


^y. 


± ^ - J- 

3 dx  ^ 2 


+ 


2 - 

y w X 

'o 


with  f(o)  = 0 and  f'(o)  = -7j/2.  This  equation  may  be  expressed  as  two  first-order 
differential  equations, 


dx 


X 


dx 


" >0 


2 _ 
osx 


wdth  f(o)  = 0 and  h(o)  = - 17.  In  this  notation 

- .2- -1/2 

(X)  = [1  - (y^UJX)  ] 

and  the  problem  is  specified  by  three  input  variables,  y^,  ~ and  t).  The  physical 
range  of  x is  0 < x < x^,  where  x^  = (y^  w)’^  is  the  point  where  y^  (x)  becomes 
infinite. 

We  can  reduce  the  problem  to  a two-variable  parameter  space  with  the 
following  transformation: 

X = y^  ojx 
f = f/y^ 
h = -h/rj 
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with  which  we  obtain 


y^(v\)  . 1 - 

— — n “ “ I 


- “ “ , - 2 
dx  (Tj/y^j)  (Vq  w) 


(Vy^)  (y^  w) 


with  f (o)  = 0 and  h (o)  = 1.  Here,  (x)  = (1  - x ) and  0 < x < 1 is  the  physical 
range  of  x.  The  problem  now  has  only  two  input  parameters,  y^  u>  and  Tj/y^’  These 
equations  are  solved  by  a standard  Runge-Kutta  integration  procedure. 

The  two  input  parameters  represent  a two-dimensional  parameter  space  which 
we  may  depict  as  a plane.  Our  first  interest  is  to  define  the  subset  of  this  plane  in 
which  self-consistent  equilibria  exist.  This  region  of  the  plane  will  be  referred  to 
as  the  equilibrium  parameter  space. 

The  density  is  given  by 


n(x)/n  = - 


’'o‘^  1 d 

^ X f (X) 
TKy^  X dx 


= y^^  (X)  h (X) 

= y^^  (X)  g (x)/g^ 


g(x)/g^  = h (x). 


The  condition  for  equilibrium  is 


-1 

g(X)  > (a; 


g(x) 


< 


-1  - 
y (X) 


for 


X ''  a 


with  the  beam  radius,  x = a,  given  by  g (a)  = y (a). 

Since  g(o)/g  = y ^ (o)  = 1,  we  can  satisfy  the  first  condition  only  if  g^  > 1. 
Furthermore,  since  y ^ (1)  = 0,  the  second  condition  requires  that  h (x)  pass  through 
zero  somewhere  in  the  interval,  0 < x <1.  Figure  7-21  illustrates  the  two  types  of 
equilibrium  solutions  which  are  possible.  Cur\'e  "a"  has  the  properly  h (x)  < y^  (x) 
for  all  X.  Consequently,  any  x in  the  range  0 < x < "’hero  h = 0,  can 

be  the  self-consistent  beam  radius  witli  tlie  suitable  choice  of  g^, 


g = [h  (X)  y (X)  ] 
o i. 

Curve  "b",  on  the  other  hand,  allows  self-consistent  solutions  with  beam  radii  in 

the  range,  x . < x < x , where  ii  (x  ) = y ^ (x  • ) and  h (x  ) = 0.  Any  x 

min  max  mm  j.  min  max 

in  this  range  can  be  the  self-consistent  beam  radius  with  g^  given  above. 

For  X > a,  we  must  calculate  h (x)  from  the  vacuum  field  equations,  with  con- 
tinuity at  X - a imposed  as  a boundary  condition  on  the  potentials  and  their  first  deriva- 
tives. The  resulting  expression  is 


h (X  > a)  = h (a)  - 


2 (y^u:)(b/y^) 


1 - a f (a)l 


_2 

a 


- £n  (^) 


which  must  pass  through  zero  somewhere  in  the  internal,  a < x < 1,  as  noted  above. 


If  we  fix  y u and  vary  rj/y  . we  always  find  equilibrium  solutions  for  small 
values  of  rj/y  • As  we  increase  rj/y  , there  is  invariably  some  value  at  uhich  the 
solution  no  longer  satisfies  h (x)  = 0 somewhere  in  the  0 < x < 1 interv'aJ.  Equili- 
brium solutions  exist,  therefore,  onh’  for  less  than  a maximum  value  which  de- 

pends on  y u-  The  region  of  the  parameter  space  where  equilibria  can  be  found  is 
shown  in  Figure  7-22.  For  y a,>  < 0.3  the  equilibrium  parameter  space  is  bounded 

bv  the  cur\’e  p/v  - 2 y a; , which  is  just  the  familiar  result, 

' 'o  o 


For  larger  values  of  y a;,  the  boundary  curve  peaks  at  r]^y^  = 0*63  and  “ 0.45, 

and  then  falls  off,  approaching  the  curve  ~ Tq  W > 2.  In  this  region 

the  maximum  value  of  rj  is  independent  of 

The  equilibrium  parameter  space  shown  in  Figure  7-22  includes  all  possible 

equilibria.  For  application  to  an  energy  storage  device,  those  equilibria  which  corres- 
pond to  high  u/y  beams  with  large  beam-field  energy  compared  to  the  guide  field  energy 
are  of  particular  interest.  It  is  necessary',  therefore,  to  examine  the  characteristics 
of  the  equilibria  which  are  possible  in  various  regions  of  the  equilibrium  parameter 


space. 


Having  fixed  y a!  Tj/y^  within  the  equilibrium  parameter  space,  we  can 


solve  the  field  equations  for  f (x)  and  h (x),  and  we  then  know 


n W.  = 


y^  (X)  h (X) 


= f (X) 


= 1 - X f (X) 


At  any  x we  may  calculate  - [h  (x)  (x)  1 and  if  > 1.  ^e  may  designate 

the  chosen  x as  the  (normalized)  equilibrium  beam  radius,  a,  which  is  numerically 

equal  to  v v (a)/c.  For  x > a , the  fields  are  given  by 
0 0 

E (X)  = (a)  = = ^ (X) 

o 

B (X)  = B 
z o 

where  B /B  = 1 - a f (a)  is  the  ratio  of  the  external  guide  field  to  that  at  the  center 
o i 

of  the  beam.  This  ratio  measures  the  diamagnetic  strength,  1 - 

In  terms  of  these  quantities,  we  may  calculate  v/ for  the  beam, 

2 a 

-B-  = — f 2Trrdrn(r) 

V ^ / 

o V me  •'o 
o 


2 (y^  w) 

If  the  beam  is  enclosed  in  a concentric  conducting  shell  of  radius,  b,  the  self-field 
energj’  per  unit  length,  W,  is  given  by 


2y  - 1 

W _ 2 'o 

M a^  (b/a)^  < V®i^  ^ 

where  W s Trb^  (B  ) is  the  guide  field  energy  per  unit  length.  The  integral  in 

MO 

this  equation  is  computed  numerically  along  with  the  field  equations  for  f and  h.  The 
potential  well  depth,  associated  with  the  beam  is  given  by 


/ 


a _2  _2_2_ 

xdx  f (X)  + a f (a)  £n  - 


07 


1 


$ 

o 

y B b 
o o 


a 


(b/a)  (B  /B  ) 
o i 


f (X)  dx  + a f (a)  fn 


b 

a 


and  this  integral  is  also  computed  numerically. 

To  complete  the  specification  of  a particular  equilibrium,  we  need  to  specify 
7^  and  either  n^  or  B^.  Knowing  y^,  we  determine  n and  hence  tc  'V  , 


cc 


D 


B 

1 


n 

o 


B 

o 


which  relates  n^  and  B^.  Since  - ( u'pL^/c)-,  we  can  determine  L^,  and  all  equili- 
brium quantities  are  then  knovm  explicitly. 

To  illustrate  the  above  discussion,  we  construct  a particular  equilibrium  in 
detail.  For  example,  for  y^  = 0.3  and  rj/y^  = 0.5,  w'hich  is  within  the  equilibrium 
parameter  space,  we  can  discuss  the  self-consistent  equilibrium  having  normalized 
beam  radius,  a = y^  v^  (a)/c  =0.8.  In  this  example  the  results  of  the  numerical  in- 
tegration of  the  field  equations  are  f (0. 8)  = -0.  826,  h (0.  8)  = 0.  22S, 

and  $^  = 0.2875  y^B^b,  for  b/a-=  2.  The  appropriate  g 
for  this  beam  radius  is  g^  = 2.63  since  y^  (0.8)  = 1.67.  With  these  results  W'e  can 

derive  t//y  = 5. 80,  u)  a/c  = 4.33,  ^ /y  =0.185  and  B/B  =1.66. 

P poo  o i 

For  definiteness,  we  now  specify  B = 20  kG  and  v =11.  Then 

10-1  13  ° -3  ° 

tOp  - 6.53  X 10  sec  , n^  = 1.5  x 10  cm  , = 0.75  cm,  and  the  beam  radius  is 

a = 2 cm  with  b = 4 cm  for  the  conducting  outer  wall.  The  guide  field  energy  is 

therefore  the  field  energj'  associated  with  the  beam  is 

W = 2. 2 kJ/cm  since  = 27.6.  The  average  energy  density  associated  with  the 

beam  field  is  44  J/cm  , and  the  depth  of  the  potential  well  between  the  bean\  center 

and  the  conducting  wall  is  75. 5 MV.  The  electrostatic  field  is  maximum  at  the  beam 

edge,  where  E = -33  MV/cm  and  falls  to  E = -16.5  MV/cm  at  the  wall 
r r 
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14  -1 

The  beam  line  density  for  this  equilibrium  is  N = 2. 3 x 10  cm  , corres- 
ponding to  36.5  fiC/cm  of  stored  charge  or  a beam  current,  1=1.1  MA.  Since 

y y = 18.3,  which  is  8.85  MeV  electron  energy,  the  total  kinetic  energy  of  the  beam 
o i 

is  320  J/cm.  The  equilibrium  fields  and  density  profile  are  plotted  in  Figure  7-23. 

On  this  figure  we  see  that  the  beam  density  is  close  to  uniform  our  to  ihe  beam  edge. 

Although  the  input  parameters  and  Ti/v^  are  useful  for  simplifying  the 

solution  to  the  field  equations,  they  are  not  the  most  useful  set  when  searching  for 
a range  of  parameters  in  which  favorable  equilibria  for  energy  storage  are  found. 

As  the  previous  example  illustrates,  there  is  no  question  about  the  existence  of  equili- 
bria which  can  store  large  quantities  of  energy.  The  remaining  question  concerns  only 
the  range  of  equilibrium  parameters  in  which  such  equilibria  exist. 

Within  the  present  formalism  we  must  specify  six  parameters,  as  V 

V V fa)/c.  V , b/a  and  either  n or  B , to  completely  specify  an  equilibrium  solution. 
e o o o 

A more  relevant  set  to  the  design  of  experiments  would  replace  the  first  two  param- 
eters, V 7)  and  Tf/y  . by  m a/c  and  a-'  h 0.  (there  are  many  other  possibilities). 

'o  o p p O 1 

The  relationship  between  these  two  sets  of  parameters  follows  from 


o y u) 
o 


or 


u a 

p 


o) 


>'o^i 


= a 


T7/y 


VoW 
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These  equations  give  us  a procedure  for  inverting  our  formalism  so  that  a/c  and 

tu  /v  f?  can  be  specified  instead  of  y uj  and  ri/y  • We  first  note  that  since  g > 1 
'"p  'o  i 'o  ' 'o  o 

for  equilibrium, 


_Q_  = 


o y n 

'o  i 


is  required.  From  Figure  7-22,  there  is  a bounded  region  of  y^  u;  for  which  this  in- 
equality is  satisfied.  We  need  only  examine  values  of  y^  us  'Within  this  range.  In  par- 
ticular, for  a)  /y  0 > 0.63,  there  are  no  equilibria, 

p o 1 _ 

Having  specified  a/c,  Wp/y^  f?.  and  a = y^  v^  (a)/c,  we  pick  a value  for 
y ol'  within  the  range  given  by  T)/y^  > Wp/Vg  Figure  7-22  and  solve  for  v^y^. 


y to  / to  a 

o / p 


From  Figure  7-22,  we  find  the  maximum  value  of  (rj/y^)  at  the  particular  choice  of 

y /jj.  For  equilibria  to  be  possible  at  all  with  these  parameters,  we  must  have 
o 


(t) 


> > 


^o^i 


If  these  conditions  are  violated,  we  choose  a new  value  for  to  and  begin  again. 

When  the  inequality  is  satisfied,  we  integrate  the  field  equations  to  x = a,  and  solve 

for  g . If  g > 1 is  satisfied,  we  have  a solution  which  has  the  correct  value  of  the 
0 0 _ 

product,  ( (Up  a/c)  Oj)-  We  repeat  this  process  for  various  values  of  y^  u-' 

until  a solution  is  obtained  which  also  has  the  correct  values  of  o)  a/c  and  a*  /y  0 

P p O 1 

individually. 


Since  Tj/y^  £ 2 tc  for  all  points  in  the  equilibrium  parameter  space,  we 


know  that  equilibria  are  possible  onlj'  for 


c y f? . 
o 1 


We,  therefore,  define  the  composite  parameter 


( a'  a/c)  (w  /y  O.) 


2y^Vg(a)/c 


In  terms  of  this  parameter,  the  two-mass  approximation  yields  the  beam-associated 
field  energy  per  unit  length  as 


— = 2 (2y  ^ - 1)  (g) 
W o \b/ 

M 


(>  V /c) f 
O 9 


, 2 


1 . (y^V^/c)  ^ 


1 b 

- + lT\  - 
4 a 


In  Figure  7-24  we  plot  this  expression  for  comparison  with  the  computed  material 


results  for  W/\\'  ,,  with  b/a  = 2.  For  y (a)/c  =0.2  and  0.5,  we  see  that  the 
M o 0 


numerical  results  approach  the  two-mass  result  as  ^ nears  unity.  This  behavior  is 


2 2 

not  surprising  since  £ = 1 corresponds  to  the  familiar  relation,  w = u!  /2  y fl  . 

^ p o 1 


As  we  increase  y v (a)/c  to  0.9,  where  the  two-mass  approximation  is  not  expected 
o 0 


to  hold,  we  find  that  W/W,,  becomes  larger  than  the  two -mass  result. 

M 


Similarly,  Figure  7-25  presents  the  numerical  results  for  i>/y  , which  in  the 

2 ° 

two-mass  approximation  is  simply  v/y  = 0.25  (^;  a/c)  . Again,  we  find  that  for 

o p 


y (a)/c  = 0. 9,  the  numerically  calculated  value  of  v/y  lies  above  the  two-mass 
o 0 o 


result. 


These  calculations  have  shown  that  there  is  no  new  limit  to  the  stored  field 


energy  as  y v /c  nears  unity.  In  fact,  self-consistent,  fully  relativistic  equilibria 
o 0 


at  large  values  of  y v (a)/c  store  larger  quantities  of  energy  than  the  two-mass 

O H 


approximation  would  predict. 
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Figure  7-24.  Beam-Associated  Field  Energy  Normalized  to  Guide  Field  Energy. 

Rostoker-Hieronymus  (Two-Mass)  Result  Plotted  for  Comparison. 
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7.1.4  Relativistic  Fluid  Equations  for  Non-Neutral  Plasma 

A majoi-  obotaclc  in  obtaining  theoretical  results  for  relativistic  non-neutra! 
systems  has  been  the  nonexistence  of  a set  of  fluid  equations  with  a physical  basis 
for  closing  the  moment  hierarchy.  In  the  absence  of  such  a set  of  equations,  the 
only  available  tool  for  theoretical  equilibrium  and  stability  analysis  has  been  the 
coupled  Maxwell-Vlasov  system  of  equations.  This  system  yields  analytical  solu- 
tions only  in  the  simplest  geometries.  Fully-relativistic  cylindrical  rigid-rotor 
equilibria,  for  example,  require  numerical  integration  for  the  equilibrium  fields. 
Stability  analysis  requires  that  the  equilibrium  particle  orbits  be  employed  as 
characteristics  to  determine  the  perturbed  distribution  function,  which  in  turn  must 
be  Integrable  over  momentum  space  to  obtain  a dispersion  relation.  Analytical 
solutions  are  therefore  very  limited  since  the  equilibrium  fields  themselves  are 
often  not  known  analytically,  and  even  then  the  particle  orbits  are  usually  known 
only  as  quadratures,  a form  which  is  not  useful  as  a characteristic.  Numerical 
solutions  of  the  Maxwell-Vlasov  system  are  possible,  but  generally  require  difficult 
and  expensive  particle  simulation  techniques,  which  are  beyond  the  scope  of  all 
but  the  largest  computers. 

By  contrast,  non -relativistic  neutral  plasma  physics  has  long  benefited  from 
the  existence  of  a closed  hierarchy  of  MHD  fluid  equations,  together  with  a number 
of  adiabatic  laws  and  an  energy  principle  for  studying  stability.  This  impressive 
theoretical  formalism  has  allowed  the  study  of  neutral  plasma  in  such  complicated 
configurations  as  symmetric  and  asymmetric  toroids,  including  toroids  of  non- 
circular cross-section.  Also,  much  of  our  understanding  of  the  physical  mechanisms 
which  underlies  gross  neutral  plasma  behavior  comes  from  fluid  analyses. 

We  have,  therefore,  undertaken  the  construction  of  a closed  set  of  fluid 
equations  for  non-neutral  relativistic  plasma.  Aside  from  the  obvious  application 
of  such  a set  of  equations  to  the  electrostatic  torus,  a fully-relativistic  system 


of  fluid  equations  would  provide  a significant  tool  in  studying  the  physics  of  diodes, 
where  the  electron  beam  is  usually  non-neutral. 

The  most  logical  place  to  begin  such  an  undertaking  is  with  the  two-mass 
approximation.  Here,  relativitj'  enters  in  only  a peripheral  way,  and  problems 
associated  with  the  non -neutrality  of  the  system  can  be  isolated  and  studied.  To 
display  the  full  complexity  of  the  general  problem,  however,  we  w'lll  discuss 
the  moments  of  the  fully-relativistic  V'lasov  equation  before  specializing  to  the 
two-mass  approximation,  Wliile  it  is  Inadequate  for  diode  problems,  the  two-mass 
approximation  is  quite  sufficient  for  describing  the  electron  beam  in  an  electro- 
static torus,  and  for  studying  relativistic  injection  techniques. 

The  fully  relativistic  Vlasov  equation  may  be  written  as 
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where  £ = ymy  and  y = [l  + p /m  c ] = [l  - v /c  j .To  compute  the 

moments  of  this  equation  we  define 
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The  relativistic  moment  equations  are  then  given  by 
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The  equations  are  the  usual  continuity  equation,  which  is  the  zeroth-moment 


of  Vlasov's  equation,  together  with  the  v-moment 


.y^Ex  [ 


Vlasov  Equ. 


, which  is 


necessary  to  determine  <v>  in  the  continuity  equation.  The  remaining  equations, 

-1 


which  are  the  moments  with  respect  to  y , v/y,  v v and  v v/y,  are  all  required 
to  close  the  second  (v-moment)  equation.  The  moment  hierarchy  which  is  obtained 
by  this  procedure  is  seen  to  be  a multiply -branched  system,  in  which  each  moment 
couples  not  only  to  the  next  moment,  but  to  several  higher  moments  as  well. 

This  property  precludes  closure  of  the  set  in  the  usual  manner,  1.  e.  by  neglecting 
a single  quantity  such  as  heat  flow.  The  situation  is  more  complicated  here  compared 
with  the  non-relatlvistic  moment  equations  because  of  the  occurence  of  y in 


relativistic  mechanics.  In  the  non-relatlvistic  domain  only  moments  with  respect 


to  1,  V,  V V,  etc.  need  be  computed  because  £ = mv  and  kinetic  energy  is  E = mv  /2, 


1.  e. , half  the  trace  of  the  v v -moment.  Relativistic  mechanics  has  £ » ym v 
2 


and  E =yTnc  , which  requires  a distinction  between  £ and  v moments,  as  well  as 
requiring  separate  moment  equations  for  energy  (the  y-moment)  and  pressure 
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(the  £ v-moment).  At  the  present  time  there  is  no  general  physical  technique  for 
terminating  the  relativistic  hierarchy  of  moment  equations. 


For  various  specific  applications  the  hierarchy  can  be  closed  by  postulating 

the  for'ni  of  thp  Histrihirtinn  function,  which  is  then  employed  to  explicitly  calculate 

the  form  of  the  unknown  quantities  required  for  closure.  This  technique  was  employed 
(2) 

by  Toepfer  , who  used  a relativistic  Maxwellian  to  study  steady-state  solutions,  and 
(3) 

by  Mosher'  , who  employed  a mono -energetic  distribution  to  study  e-beam  interac- 
tions with  high-Z  plasma.  The  results  obtained  by  this  technique  are  only  as  valid  as 
the  assumed  distribution  function.  In  Toepfer's  work,  for  example,  the  use  of  a 
Ma.x\\’ellian  to  study  steady-state  problems  is  reasonable  since  the  departure  from  a 
Maxwellian  behavior  is  only  valid  in  the  collislonal  regime  where  equilibrium  times 
are  short. 

We  turn  now  to  the  two-mass  approximation,  in  which  the  electrons  are  rela- 


tivistic in  one  direction  only,  i,  e. , v = V + (v  - V),  where  V » v , v - V . In 

^ ^ 2 -1/2'  ^ 

this  approximation,  V is  a constant,  8^sV/c,  and  = (1  “ ) • Rigid-rotor 

Vlasov  equilibrium  and  stability  analysis  in  both  cylindrical  and  large  aspect-ratio 

toroidal  geometries  have  been  performed  in  this  approximation  by  Rostoker  and 
141 

Hieronymus  . These  results  have  been  presented  in  a previous  Final  Report. 

The  Vlasov  equation  in  the  two-mass  approximation  is  written  as 
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where  f = f(x,  v,  t)  and  i , n refer  to  directions  perpendicular  and  parallel  to  V. 


A.  J.  Toepfer,  Phys.  Rev.  1444  (1971). 

D.  Mosher,  Bull.  Am.  Phys.  So.  1^,  856  (1974). 

J.  Hieronymus,  Ph.  D.  Thesis,  Cornell  University  (1971). 


We  define  the  quantities. 
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With  these  definitions,  the  first  three  moments  of  Vlasov's  equation  are 


9n 


— + V • n <v>  = 0 

ot 


d<v> 

nm  — - — = -V 


dt 


^ ne  r 1 


ne 

3 


E + - (<v>  X B) 
-ii  c - -1 


T T 

> 


^ P + V • |^< V > P + qJ  + P • V < V>  + I^P  . V < V>J 

i]  + [l  ll  " ° 

o y me  ^ J 

o 


where  d/dt  = S/3t  + <v>  • V and  the  cross-product  operators  (Bx)i,  and  (xB)i  are 
defined  by 


w (xB)i,„  = (wtB)x,„ 

for  any  vector  w.  Here  again  a,  h are  with  respect  to  V.  The  quantity  [P  • 7 <v>]'^ 
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is  the  transpose  of  P • V <v>. 

This  set  of  equations  may  be  closed  by  setting  the  heat  flow  tensor  to  zero, 
i.e.  Q = 0.  Unlike  ordinary  one-fluid  MHD  theor>’,  however,  it  is  unnecessary" 
to  assume  a form  of  Ohm's  law  to  complete  the  equations.  In  particular,  the  infinite 
conductivity'  assumption,  E + <v>  x B/c  = 0,  is  very  restrictive  for  the  single- 
species  equations  considered  here.  This  assumption  would  drop  all  electromagnetic 
terms  from  the  force  equation  since  for  a pure  electron  plasma  the  current  density 
is  J = -ne  <v>  and  the  charge  density  is  p = -ne.  The  infinite  conductivity  equation 

( 5 

is  therefore  equivalent  to  setting  dE  + J x B/c  = 0.  In  the  one-fluid  MHD  equations 
the  situation  is  entirely  different  because  infinite  conductivity  is  expressed  as 
E+  U X b/c  = 0 where  _U  is  the  center-of-mass  velocity  given  by 

_ mn  <y> 

— r mn 


and  the  summation  is  over  particle  species. 

To  check  the  fluid  equations,  we  begin  by  comparing  the  steady -state  solutions 
from  the  Vlasov  formalism  with  the  steady-state  fluid  equations.  To  construct 
an  exact  tw'o-mass  Vlasov  equilibrium,  we  begin  with  the  two-mass  equations  of 
motion. 
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See  for  example  l.B.  Bernstein  and  S.K.  Trehan,  Nucl.  Fusion  1,  3 (1960). 
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for  a cylindrical  beam  with  self-consistent  fields  E^,  which  are  functions 

of  the  radial  coordinate  only.  The  constants  of  motion  associated  with  the  equations 
of  motion  are 
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We  write  the  equilibrium  distribution  function  as 
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The  moments  of  F(x,  v)  which  are  required  to  determine  the  equilibrium 
fields  are 
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where  © (x)  Is  the  unit  step  function,  which  is  unity  for  x > 0 and  zero  for  x < 0. 

The  beam  radius,  r = a,  is  given  by  ji(a)  = 0,  where 

•i 

, , 1 3 e 2 

l'(r)  = -y  in  (V  + — A ) + U(r) 

<2  0 o Z 

y me 
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and  we  have  assumed  ib{r<a.)  < 0 and  i/)(r>a)  > 0 as  can  be  demonstrated  a posteriori . 

Higher  moments  of  F (x,  v),  which  are  required  in  the  steady-state  fluid  equations  ^ 

are 

^r  = ^ee  = (a-r) 

Q = 0 . 

The  remaining  components  of  the  pressure  tensor  are  zero,  and  also  <v  > = 0. 

r 

We  can  now  show  that  the  steady-state  fluid  equations  are  satisfied.  The 
continuity  equation  is 

i 

V • n < v>  = n V . < v>  + <v>  • vn  = 0 | 

I 

which  is  satisfied  because  <v^>  = 0 and  all  quantities  depend  only  on  r.  We  therefore 
have 

n V . <v>  = <v>  • vn  = 0 . 

\, 

The  fl  and  z components  of  the  momentum  equation  are  trivially  staisfied.  The  1 


li  radial  component  is  also  satisfied  because 


The  pressure  equation  is  satisfied  because 
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We  used  V • <v>  - 0,  <v^>  - 0 and  P = P (r)  in  obtaining  the  above  results.  The 
steady-state  fluid  equations  are  therefore  satisfied  by  the  moments  of  a Vlasov 
equilibrium. 


To  complete  the  equilibrium  calculation,  we  compute  the  fields  from  Maxwell's 
equations.  The  result  can  be  used  to  calculate  <v  > explicitly, 


<v  > = V + 
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A (r)  = V I 
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where  I is  the  modified  Bessel  function  of  zero  order  and  cc  = 47m  e /v  m. 

I P I o ' 

For  this  solution  to  be  meaningful,  we  must  require  V » <v  > - V , which  implies 
that  the  argument  of  the  Bessel  function  be  small.  The  equilibrium  is  therefore 
valid  when 

3 

I'/y  « 1/4 


which  is  only  a weak  restriction  on  v/y^.  In  this  limit  the  equilibrium  obtained 
above  reduces  exactly  to  the  Rostoker-IIieronymus  two-mass  equilibrium. 
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Returning  now  to  the  full  set  of  fluid  equations,  closed  with  the  assumption 
2 = 0,  we  can  establish  the  adiabatic  law  for  isotropic  pressure,  P = p i , by 
taking  the  trace  of  the  pressure  equation, 

3 if-  + 5o  V . <v>  = n , 

UL 

We  use  V • <v>  = -n  ^ dn/dt,  which  follows  from  the  continuity  equation,  to  obtain 
the  adiabatic  equation  of  state 


_d 

dt 
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0, 


The  same  result  follows  from  ordinary  one-fluid  MHD,  provided  the  additional 
assumption  of  infinite  conductivity  is  made.  The  adiabatic  law  for  the  two-mass 
fluid  equations  does  not  require  this  assumption.  As  an  aside,  we  note  that  the 
double -adiabatic  law-,  or  Chew-Goldberger-Low  theory,  does  not  follow  from 
our  set  of  fluid  equations. 

Using  the  adiabatic  law  derived  above,  the  linearized  fluid  equations  may 
be  written  as 
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where  we  have  dropped  the  angle  brackets,  \ from  our  notation.  The  presence 

of  the  zero-order  velocity,  v , is  required  by  Eamshaw's  theorem,  which  states 

~o 

that  there  are  no  static  equilibria  for  an  unneutralized  collection  of  charged  particles. 
A non-neutral  equilibiium,  therefore,  must  be  dynamic  in  character.  This  feature 
of  the  equations  spoils  the  usual  deri ration  of  the  energy  principle,  which  does  not 
apply  in  its  conventional  form  to  non-neutral  systems. 

To  study  the  stabiliU’  of  finite -radius  equilibria,  such  as  the  rigid -rotor 
equilibrium  described  earlier,  the  linearized  fluid  equations  become  a coupled  set 
of  ordinary  first-order  differential  equations  which  in  general  may  be  solved  numeri- 
cally by  the  Runge-Kutta  integration  procedure.  Anahlical  solutions  are  possible 
in  special  cases. 

As  an  example  we  analyze  the  well-known  diocotron  instability',  and  compare 

our  result  with  Levy's  non-relativistic  analysis.  We  assume  an  equilibrium  specified 

o o o o o 

by  n (r)  v„  (r),  v (r),  p (r),  E (r),  B (r),  B (r)  for  an  annular  beam  with 
o z o r p z 

inner  radius,  a,  and  outer  radius,  b.  The  diocotron  mode  is  an  electrostatic, 
incompressible,  azimuthal  mode  of  an  unneutralized  beam.  We  therefore  specialize 
the  analysis  to  perturbations  which  satisfy 
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and  assome  all  perturbed  quantities  to  be  of  the  form 
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with  no  dependence  on  the  2 -coordinate. 
We  define  the  frequencies 
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and  solve  for  v in  the  linearized  momentum  equations. 
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Using  the  linearized  pressure  equation,  we  obtain 
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The  perturbed  density  then  follows  from  the  linearized  continuity  equation 
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which  is  the  source  function  in  Poisson's  equation. 
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For  an  equilibrium  with  constant  v and  uniform  density,  n , within  the  annulus  and 

z o 


for  a pressure  profile,  P^lr),  is  goes  smoothly  to  zero  at  the  beam  edges,  Poisson's 


equation  is  simply 


t 
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7.1.5  Switching  an  Unneutralized  Electron  Beam  in  a Torus 


Any  analysis  of  switching  an  unneutralized  electron  beam  in  a torus  must  assume 
linear  stability  of  the  beam.  Without  linear  stability  the  energy-  storage  would  not  be 
possible  since  nearly  any  pure  electron  mode,  if  unstable,  would  grow  at  a rate  com- 
parable to  ajpp  = (477  ne^/y^m)^'^^.  At  n = 10^^  cm”^  and  y^  = 5,  we  have  — 

4 X 10”^^  sec.  Any  beam  instability  would  e-fold  in  amplitude  in  <1  ns  at  this  density 
and  energy.  The  rapid  growth  rate  for  instabilities  suggests  one  approach  to  switching, 
viz.  changing  some  parameter  to  trigger  an  instability. 

Given  linear  stability,  the  only  instabilities  possible  are  those  which  result 

from  a large-amplitude  perturbation,  one  comparable  in  magnitude  to  the  equilibrium 

value  for  some  quantity.  For  example,  by  inserting  a rod  into  the  torus,  we  can 

produce  a kink  perturbation  with  initial  amplitude  comparable  to  the  beam  radius. 

Alternatively,  Rosenbluth  has  suggested  that  we  switch  the  torus  by  injecting  a fast 

plasmoid  into  the  beam,  thereby  rapidly  charge  neuti-alizing  a local  region,  i.e, , 

6E  ~E  , which  results  in  an  unstable  neutral  plasma  mode, 
r r 

A separate  line  of  reasoning  takes  advantage  of  the  linear  stability  of  the  beam. 
Such  approaches  as  dropping  a pellet  into  the  beam  or  magnetically  moving  the  beam 
into  a metal  target  have  been  considered  and  rely  on  beam  stability  during  switching. 

The  basic  objective  in  all  switching  schemes  is  to  charge  neutralize  and  stop 
the  beam  in  a time  of  order  an  electron  transit  time  ~27rR/c.  The  external  circuitiy 
has  time  constants  which  are  long  compared  with  the  switching  time,  and  therefore 
the  field  energy  will  be  dissipated  in  the  target  material  heating  it  and  producing  x-rays 
by  Bremsstrahlung. 


Plasma  Switch  (Rosenbluth) 


Fast  theta-pinch  hydrogen  plasma  guns  can  be  made  to  produce  3x10 


electron-ion  pairs.  Typically,  about  1%  of  these  particles  are  emitted  in  a fast 
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plasmoid,  with  velocities  of  10  - 10  cm/sec  (plasmoid  energy  <100  J).  Since  the 
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total  number  of  electrons  in  the  beam  is  ~10^^  in  a megajoule  system  (and  in 

the  STP),  the  fast  plasmoid  will  charge  neutralize  a local  region  of  the  beam  in 
~ 10  ns. 

A neutral  plasma  confined  in  a torus  cannot  carr>’  current  in  excess  of  the 

Kruskal-Shafranov  limit.  A fully  neutralized  electrostatic  torus  would  therefoi-e  be 

grossly  unstable.  With  the  disappearance  of  the  electrostatic  field  in  a local  region, 

the  B field  will  be  unopposed  and  cause  the  beam  to  pinch.  This  unstable  pinch 
6 

motion  can  stop  a high  v/y  beam.  Targets  on  the  wall  can  be  used  for  conversion  to 
x-rays,  or  alternatively,  the  plasma  itseK,  which  is  strongly  heated  during  the  beam 
pinching,  can  become  the  x-ray  source. 

"Rod-Plunger"  Switch 

By  making  a perturbation  in  the  wall  geometn-,  we  affect  the  boundan  conditions 
which  determine  the  equilibrium  fields.  It  is  possible  to  actually  destroy  the  equilibrium 
locally  in  this  manner  and  force  a specific  perturbation  (i.e.  a kink)  to  grow.  We 
imagine  a cylindrical  rod  which  is  inserted  through  the  torus  wall  and  pushed  toward 
the  minor  axis.  Initially,  the  rod  simply  compresses  all  the  fields,  E , B^,  B about 

equally.  When  the  rod  has  been  inserted  to  a depth  comparable  to  its  diameter, 

however,  the  magnetic  field  lines  will  begin  to  slip  around  the  sides  of  the  rod,  while 

the  electric  field  lines  must  continue  to  terminate  on  the  rod.  Further  insertion  of 

the  rod,  therefore,  results  in  an  enhancement  in  E^  , but  little  further  compression 

of  B„  or  B . 

B z 

Increasing  E causes  the  azimuthal  drift  speed  v’  to  increase  as  electrons 
r 6 

pass  near  the  rod,  forcing  a kink-shaped  perturbation  to  form.  Let  6E^  be  the  change 

in  E due  to  the  rod.  Then  the  azimuthal  speed  is 
r 


V^B 
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impossible  and  the  electrons  will  be  freely  accelerated  toward  the  rod. 


The  switching  process  can  be  triggered  by  biasing  the  rod  positive  as 

nears  unity.  A bias  voltage  -10-100  kV  can  be  rapidly  applied,  and  results  in  the 

iormation  of  a large-amplitude  perturbation  (amplitude  —beam  radius)  in  a time  1 ns. 

The  large  kinl-;  perturbation  then  grows  very  rapidly,  with  a characteristic  grov,ih 

time  <1  ns,  and  drives  the  beam  into  the  rod  which  acts  as  a target. 

This  mechanism  does  not  rely  on  rapid  insertion  of  the  rod,  but  only  on  fast 

onset  of  the  instability.  The  large  E and  B fields  of  the  beam  itself  are  utilized  in 

r 0 

switching,  the  bias  voltage  acting  only  as  a trigger. 

Pellet  Switch 

Perhaps  the  simplest  approach  to  switching  consists  of  dropping  a metal  pellet 

into  the  beam  to  charge  neutralize  and  stop  the  beam.  The  pellet  then  heats  up  and 

emits  X rays.  The  main  problem  is  getting  the  pellet  into  the  beam.  By  charging 
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the  pellet  positive,  the  large  potential  between  the  wall  and  beam  surface  (-10  - 10  \ ) 
can  be  used  as  an  accelerator.  Pellet  velocities  up  to  -10  cm/sec  are  attainable  in 

—8 

this  way.  While  the  total  beam  mass  (-10  g)  is  very  small,  the  B^  and  B^  fields 
prevent  the  beam  electrons  from  being  drawn  into  the  pellet.  Clearly  a technique  for 
moving  the  beam  into  the  pellet  is  required. 

A means  of  shifting  the  beam  equilibrium  through  the  application  of  a vertical 
magnetic  field  appears  feasible.  This  approach  is  discussed  in  detail  in  the  following 
section.  A vertical  magnetic  field  of  ~300  G suffices  to  displace  the  beam  a distance 
comparable  to  its  diameter  in  -10  ns.  A description  of  a technique  for  producing  such 
a field  in  -10  ns  is  also  presented. 

Switching  By  Vertical  Magnetic  Field 

A vertical  magnetic  field  is  generally  employed  to  adjust  the  equilibrium 
position  in  a Tokamak.  We  consider  this  possibility  for  a high  v/y  electron  Ix'am  in 
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a torus.  The  objective  is  to  displace  the  electron  beam  a distance  of  the  order  of 
the  diameter  of  the  beam  in  a time  of  the  order  of  10  ns.  We  conclude  that  this  is 
feasible  with  state-of-the-art  electrical  engineering. 

We  begin  with  the  orbit  equations  for  relativistic  electrons  in  a large  aspect 
(11 

ratio  torus 
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(x,y,z)  are  local  cartesian  coordinates  as  illustrated  in  Figure  7-26,  (r,  6)  are  local 

polar  coordinates  corresponding  to  (x,y).  These  equations  are  identical  to  Eq.  (11) 

of  reference  1 except  for  the  addition  of  an  external  magnetic  field  B^  which  results 

in  the  additional  term(e/c)v  B in  the  x-component  of  the  orbit  equation.  Since 
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-e(E^  - f 


my  X 
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where  = 47rne^/y  m is  the  beam  plasma  frequency,  the  toroidal  term  my  V“/K 
p o 

can  be  removed  by  the  transformation 


Rostoker,  Particle  Accelerators  5,  93  (1973). 


This  means  that  the  beam  is  not  centered  in  the  torus  but  shifted  to  the  position  B a 
2 2 2 

distance  2y^  \ ^ center.  Similarly,  it  is  evident  that  the  effect  of 

the  vertical  field  is  to  shift  the  center  of  the  beam  from  B to  B'  by  a displacement 
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where  n = eB  /y  me.  For  example,  if  y =21,  n = 10  cm  and  B = 300G, 
y y o ^o  ’ y ’ 

Ax  = 4.4  cm.  For  the  STP,  we  have  n = 4 x 10^^  cm“^  and  y = 5,  which  yields 


B = 195  G for  Ax  = 4 cm. 

y 

The  X and  y components  of  the  orbit  equations  can  be  written  as  follows: 


X = f2  x-n  y + V n (t) 
z z y 


y = n y + n X 
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where 
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In  terms  of  the  complex  variable  = x + iy  these  equations  can  lx?  combine 
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Ct^'»  C<0)  initial  values  for  ^(t)  and  o;^  =f)  /f)^  is  the  slow  precession  frequency. 

It  has  been  assumed  that  to  « 0 or  0^/0?  «1.  The  orbital  motion  descriljed  by 

s z ^ 

equation  is  illustrated  in  Figure  7-27. 

The  term  e*^z^  describes  the  fast  precession  at  the  gyrofrequency  of  the  large 
toroidal  field.  The  gj'ro radius  of  the  precession  is  quite  small.  It  is  of  order 
(velocity /f)  ) <c/f)  . The  term  exp  i to  t describes  a slow  precession  about  a center 

z z 2 s 

that  is  shifted  by  V f)  /O  compared  to  the  case  Q =0.  This  agrees  with  our 

2 y y 

expression  for  Ax  above.  However,  it  is  now  clear  that  the  timescale  for  establish- 
ing this  deflection  is  determined  by  the  precession  frequency  to  ; the  speed  with  which 

s 

the  beam  moves  from  B to  B'  in  Figure  7-26. 
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We  assume  that  B = 300  G,  B =20  kG  so  that  Ax/At  = 4.5  x 10  cm/sec  and  the 

y z 

time  for  one  beam  to  displace  from  B to  B'  is 


At  ^ 10  ns 


which  is  also  the  result  for  STP  where  B =15  kG. 
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We  have  assumed  that  is  turned  on  instantaneously.  In  order  to  displace 
the  beam  in  about  10  ns,  it  will  be  necessary  to  increase  from  zero  to  300  G in 
about  this  time. 

We  assume  a set  of  toroidal  windings  to  produce  B^  as  illustrated  in  Figure  7-28. 
In  order  to  estimate  the  inductance  of  the  coil  we  calculate  the  inductance/unit  length 
neglecting  the  toroidal  curvature.  For  w » b 


£ — 8 log  (— j nanohenries/centimeter  . 


The  total  inductance  is 


L ^ 2nR£  ^16nR  log  (wA»)  nH 


where  R is  the  major  radius  of  the  torus.  For  R = 50  cm,  w = 30  cm,  b = 1.5  cm, 
L = 7500  nH. 

The  magnetic  field  in  the  vicinity  of  the  electron  beam  will  be 


B . 
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For  B = 300  G,  1 = 11.25  kA,  while  in  the  STP,  1 = 7.3  kA  will  be  required. 

y 

We  assume  that  the  coil  will  be  driven  by  a high  voltage  pulse  line.  In  order 
to  reduce  the  effective  inductance  the  coil  would  be  split  into  N coils  as  illustrated 
in  Figure  7-29.  Each  of  the  coils  is  to  be  connected  iii  parallel  and  driven  by  the  pulse 
line.  The  pulse  line  must  supply  a current  N1  and  the  magnetic  energy  is 


1 2 1 L 2 
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Therefore  the  load  inductance  for  the  pulse  line  is 


L = L/N^  . 
E 


The  effective  inductance  for  the  pulse  line  is  = 7500/64  = 117  nH  assuming  N = 8 
as  illustrated  in  Figure  7-29.  The  equivalent  circuit  for  the  pulse  line  is  illustrated 
in  Figure  7-30.  The  coils  can  be  represented  as  lumped  parameter  circuits  rather 
than  transmission  lines  since  the  transit  time  for  electromagnetic  waves  in  each 
coil  is  (2ffR/8c)  e:  1.2  ns.  From  analysis  of  the  pulse  line  illustrated  in  Figure  7-30, 
thq.  current  in  the  inductive  load  would  be 


where 


i(t)  = -r-T  — e 
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Time  is  measured  from  when  the  switch  S is  fired;  t is  the  transit  time  for  line  Z . 

2 2 

This  equation  is  valid  for  0 <t  - t <2t  . 2t  . Assuming  L » L and  Z a:Z  .t  »1 

2 12  E s 12  2 

and  the  equation  simplifies  to 


V -(t  - t„) 
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I(t)  = ^ 1 - e ^ 
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The  risetime  of  the  current  is  L /Z  and  since  this  must  be  about  10  ns,  Z ^ 120. 

£ 2 2 

The  current  rises  to  V /Z  which  must  be  8 x 11.25  = 90  kA  so  that  the  pulse  line 

O 2 

voltage  is  about  1 MV.  In  STP  a pulse  line  voltage  of  700  kV  would  be  required. 
Clearly,  this  is  a standard  state-of-the-art  pulse  line.  We  assume  that  an  oil 
dielectric  would  be  employed  and  the  radius  of  the  outer  cylinder  would  be  about  the 


same  as  the  torus,  i.e.  R =50  cm,  since  Z = 60  log  R /R.  the  inside  radius  would 

o o i 

be  R.  = 41  cm  to  give  an  impedance  of  120.  The  energy  stored  in  the  magnetic  field 
of  the  coil  is 
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- X 7500  X 10  X (11.25  X 10  ) 


= 474  Joules  . 


The  energy  stored  in  the  pulse  line  is  (1/2)  CV  /unit  length  where 


C = 


„ ,5  2 log  R /R. 

9 X 10  o i 


^F/cm  . 
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C = 0.28  X 10  ^F/cm  and  the  energy  stored  is  1.4  J/cm.  The  pulse  line  should 
be  at  least  20  ns  in  electrical  length  which  involves  a physical  length  of  6 m and  a 
stored  energy  of  840  J. 

It  is  possible  to  place  the  B^-coil  inside  the  toroidal  metal  wall.  However,  it 

is  much  easier  to  place  it  outside  as  indicated  in  Figure  7-28.  In  this  case,  the  metal 

wall  must  have  a poloidal  slot  or  the  magnetic  field  will  be  shielded  and  will  not 

3/2  , 

penetrate.  For  penetration  of  the  electron  beam  the  inductive  skin  depth  is  (c/ccp) 

17.5  cm  so  that  this  is  not  a problem  for  large  y^. 

7. 2 STP  EXPERIMENTAL  PROGRAM 

7.2.1  Introduction 

During  the  past  fiscal  year  construction  was  completed  on  the  STP  toroidal 
energy  storage  experiment  at  Maxwell  Laboratories,  Inc. , and  experiments  were  ini- 
tiated to  study  the  injection  and  storage  of  unneutralized  electron  columns  and  beams  in 
toroidal  geometry.  During  the  course  of  the  year,  some  modifications  were  made  to 
the  machine  to  correct  magnetic  field  errors  resulting  from  nonuniformities  in  the 
stainless  steel  vacuum  vessel.  In  addition,  during  the  latter  part  of  the  year  vertical 
magnetic  field  windings  were  added  to’  provide  the  vertical  magnetic  field  necessary  for 
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controlling  the  equilibrium  position  of  the  beam.  A plan  view  of  the  machine  as 
modified  is  shown  in  Figure  7-31,  and  a cross  sectional  drawing  is  shown  in  Figure 
7-32.  The  vertical  magnetic  field  windings  are  driven  by  a bank  of  marine  batteries 
which  is  switched  by  an  SCR  circuit. 

During  the  course  of  the  year,  a 40  kV  thermionic  electron  injector  was  de- 
veloped which  injects  ~100  pcoulombs  of  charge  into  the  machine.  This  injected 

10  -3 

charge  corresponds  to  an  electron  density  of  ~10  cm  and  creates  a potential  well 

of  ~300  kV.  This  amount  of  charge  is  sufficient  to  provide  the  E x B , particle  rota- 

r (D 

tional  transform  necessary  for  confining  accelerated  particles.  Efforts  near  the  end 
of  the  year  have  focused  on  acceleration  of  the  injected  charge  into  a relativistic  beam. 
WTiile  these  efforts  have  not  yet  been  successful,  we  believe  we  have  identified  the 
problem  area  involved  and  have  outlined  an  experimental  program  which  will  lead  to 
significant  relativistic  beam  currents. 

7. 2. 2 Diagnostics 

Two  different  kinds  of  probes  have  been  used  to  measure  the  charge  injected 
into  the  STP  machine.  Figure  7-33  shows  one  of  these  types  of  probes.  This  probe 
is  called  a wall  probe  and  consists  of  a flat  circular  stainless  steel  disk  mounted 
flush  with  the  stainless  steel  vacuum  vessel  wall  but  electrically  insulated  from  the 
vacuum  wall.  The  stainless  steel  disk  is  connected  to  a 50fl  cable  and  500  termina- 
tion through  a high  vacuum  electrical  feedthrough.  This  probe  provides  two  separate 
measurements  of  the  injected  charge.  First,  the  probe  measures  the  charge  induced 
on  the  wall  of  the  torus  by  the  accumulation  of  the  electron  column  inside  the  machine. 
Since  the  probe  is  connected  to  ground  through  500,  the  induced  charge  on  the  disk 
produces  a voltage  across  the  termination  when  it  flows  to  ground.  Integrating  the 
voltage  appearing  across  the  500  termination  gives  a measure  of  the  total  charge  in- 
jected into  the  torus. 
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TOROIDAL  ELECTRIC 
FIELD  WINDING  (E.) 


Figure  7-32.  Cross  Section  of  the  STP  Machine. 


In  addition,  the  wall  probes  couple  capacitively  to  the  diocotron  surface  flute 
modes  which  can  exist  in  an  unneutralized  electron  colunan.  The  frequency  of  the 
diocotron  mode  is  given  by 


As  a result,  the  electron  density  may  be  determined  by  measuring  the  diocotron  fre- 
quency if  the  magnetic  field  is  known.  Specifically, 

n = 6,9  X 10®  f^B 
a 

2 

where  the  diocotron  frequency  f^  is  in  Hz  and  the  magnetic  field  B is  in  Webers/m  . 
For  most  of  the  experiments  described,  wall  probes  were  placed  at  two  different 
toroidal  locations. 

Another  probe  used  is  shown  in  Figure  7-34.  This  probe  consists  of  a glass 
tube  with  a . 020"  diameter  tungsten  wire  protruding  from  the  end.  The  probe  which 
is  movable  radially  through  a vacuum  fitting  is  used  to  measure  the  depth  and  profile 
of  the  potential  well  created  by  the  non-neutral  charge  collection.  This  provides  a 
measure  of  the  density  in  the  electron  column.  The  probe  is  connected  through  a 
lOMfltolOkn  voltage  divider  to  an  operational  amplifier  which  drives  a terminated 
son  cable  to  an  oscilloscope.  For  most  of  the  measurements  described  in  this  report, 
the  potential  probe  was  located  toroidally  180°  from  the  injection  sector. 

Assuming  that  the  electron  density  distribution  in  the  electron  column  is 
uniform,  the  density  is  related  to  the  potential  well  depth  by 


V 

o 


= neb 
4e 
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where  V is  the  potential  at  the  center  of  the  column  and  b is  the  beam  radius.  Taking 
o 

b = 8 cm  (the  vacuum  wall  radius)  gives 

n = 3. 45  X 10^  V . 

o 

Thus,  the  potential  probe  provides  a third  method  of  measuring  the  electron  density 
inside  the  column. 

7.2.3  Low  Voltage  Injector  Studies 

The  first  injector  assembly  used  to  inject  charge  into  the  STP  machine  is  shown 
in  Figure  7-35.  This  injector  consisted  of  12  thoriated,  carburized,  .024”  diameter, 
4”  long  filament  wires  located  in  one  quadrant  of  the  minor  circumference.  Molyben- 
dum  plates  are  placed  parallel  to  the  filament  wires  and  the  filament  bias  voltage  is 
applied  across  the  plates.  This  results  in  an  electric  field  which  is  perpendicular  to 
the  toroidal  magnetic  field,  causing  the  electrons  to  acquire  an  ExB  drift  velocity  to- 
ward the  minor  axis. 

Figure  7-36  shows  oscilloscope  traces  of  machine  operation  using  this  kind  of 
injector.  For  this  shot  the  magnetic  field  had  a risetime  of  2 msec,  reaching  a maxi- 
mum (at  the  minor  axis)  of  4. 7 kG.  An  injector  bias  voltage  of  4.  0 kV  was  used,  and 
injection  was  begun  after  1.8  msec.  The  bias  voltage  was  crowbarred  after  200  ^xsec 
terminating  the  injection  phase.  The  lower  trace  in  the  upper  right  hand  photograph 
shows  the  potential  well  depth  time  history  at  a sensitivity  of  10. 7 kV/div.  As  can  be 
seen,  the  potential  well  depth  builds  to  a maximum  during  injection  and  then  begins  to 
decay  after  the  injection  process  is  terminated.  The  upper  trace  in  the  lower  left 
photograph  shows  the  wall  probe  signal  gives  the  injected  charge  (~8  ^coulombs  for 
this  shot)  which  is  in  good  agreement  wdth  the  potential  well  depth  measurement  if  an 

approximately  uniform  charge  distribution  is  assumed. 

The  bottom  trace  in  the  lower  right  photograph  shows  the  diocotron  oscilla- 
tions on  a fast  time  sweep  (1  ^ sec/division).  These  oscillations  occur  late  in  time 
when  the  potential  probe  is  located  near  the  center  of  the  electron  column.  (This 


Figure  7-35.  Linear  Electron  Injector  - End  View. 


indicates  that  the  potential  probe  perturbs  the  electron  column  sufficiently  to  cause 
the  diocotron  oscillations  to  begin).  Measurements  of  the  electron  density  from  the 
diocotron  frequency  give  values  of  electron  density  which  are  somewhat  lower  than 
the  densities  obtained  from  the  probe  measurements.  However,  this  is  probably 
reasonable  since  the  diocotron  waves  are  surface  modes  and  the  electron  density  is 
undoubtedly  somewhat  lower  near  the  edge  of  the  electron  column. 

In  Figure  7-37  we  show  probe  signals  for  a case  where  the  magnetic  field  rise 
time  was  3.5  msec,  the  peak  magnetic  field  was  5.8  kG,  and  the  bias  voltage  was 
5.  0 kV.  The  potential  well  depth  in  this  case  was  ~ 30  kV.  The  lower  two  traces  in 
Figure  7-37  show  diocotron  wave  oscillations  taken  with  two  different  wall  probes 
located  approximately  90°  apart  toroidally  around  the  machine.  The  diocotron  modes 
should  have  zero  phase  shift  parallel  to  the  magnetic  field.  As  can  be  seen,  the  two 
wall  probe  signals  are  in  phase.  These  measurements  thus  provide  verification  that 
the  observed  waves  are  indeed  the  diocotron  modes. 

In  Figure  7-38  we  show  the  variation  of  potential  well  depth  with  the  injector 
bias  voltage  for  the  low  voltage  injector.  As  can  be  seen,  the  well  depth  (and  hence 
the  charge)  was  found  to  scale  linearly  with  the  bias  voltage.  This  scaling  indicates 
that  the  injected  charge  can  be  significantly  increased  by  increasing  the  bias  voltage. 
Advantage  was  taken  of  this  scaling  later  in  the  year  to  obtain  an  increase  of  ~10  fold 
in  the  injected  charge  and  resulting  potential  well  depth. 

The  radial  potential  profile  for  the  low  voltage  linear  injector  is  shown  in 
Figure  7-39.  The  profile  can  be  reasonably  well  fit  by  a parabolic  curve  over  much 
of  the  minor  radius.  This  is  consistent  with  an  electron  column  density  which  is 
approximately  uniform  over  the  column  diameter. 

Figure  7-40  shows  the  variation  in  potential  well  depth  as  a function  of  the 
time  delay  between  the  start  of  the  magnetic  field  and  the  start  of  the  injection  phase 
As  can  be  seen,  the  peak  potential  well  depth  attained  is  maximum  late  in  time  for  both 
the  3.5  and  2 msec  magnetic  field  rise  times.  More  recent  results  with  another  in- 
jector of  only  slightly  different  geometry  do  not  show  this  same  variation,  however. 
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Figure  7-38.  Potential  Well  Depth  as  a Function  of  Injector  Bias  Voltage  for 
the  Low  Voltage,  Linear  Injector. 
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Figure  7-39.  Radial  Potential  Profile  for  the  Low  Voltage, 
Linear  Injector. 
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This  indicates  that  this  behavior  may  be  due  to  magnetic  field  errors  (due  to  mag- 
netic field  penetration  through  the  stainless  vacuum  liner)  which  were  accentuated 

by  this  injector  geometry.  j 

The  variation  of  the  potential  well  depth  with  the  peak  magnetic  field  at  the 
minor  axis  is  shown  in  Figure  7-41  for  both  early  and  late  time  injection.  The  poten-  | 

tial  well  depth  was  observed  to  be  only  a weak  function  of  the  magnetic  field  for  this  I 

injector. 

In  order  to  study  the  containment  time  of  the  electron  column,  provisions 

have  been  made  to  crowbar  the  injection  voltage  after  a time  ^t  , which  can  be  set 

cb 

by  an  external  time  delay  control.  After  the  injection  voltage  has  been  removed,  the 
electron  density  begins  to  decay,  and  this  decay  has  been  observed  for  a number  of 

■ 

different  conditions.  For  example.  Figure  7-42  shows  the  variation  of  the  potential 
decay  as  a function  of  the  peak  magnetic  field  strength.  Each  cur\’e  in  Figure  7-42 
has  been  normalized  to  its  own  peak  value  so  that  the  decay  rates  can  be  more  easily 

compared.  As  can  be  seen,  a systematic  improvement  in  containment  time  occurs  ; 

as  the  peak  magnetic  field  is  increased.  Similarly,  Figure  7-43  shows  the  decay  ^ 

rate  for  several  different  injection  times.  This  figure  shows  a systematic  improve- 
ment in  containment  time  as  the  injection  window  is  moved  to  later  times.  However, 
since  the  magnetic  field  during  the  decay  time  also  increases  as  the  injection  window 
is  moved  to  later  times,  it  is  interesting  to  compare  the  data  from  Figures  7-42  and 
7-43  on  a single  plot  in  terms  of  magnetic  field.  This  is  done  in  Figure  7-44  by 

plotting  the  decay  time  (time  for  the  potential  to  decay  to  1/e  times  the  peak  value)  i 

against  the  magnetic  field  at  the  start  of  injection.  The  upper  curve  is  obtained  from 
the  data  of  Figure  7-43  where  the  time  at  which  injection  starts  is  varied,  and  the 

lower  curve  is  obtained  from  the  data  of  Figure  7-42  where  the  peak  magnetic  field  1 

is  varied.  A primary  difference  between  the  two  curves  is  that  upper  curve  was  taken 
with  a bias  voltage  of  4. 0 kV  while  the  lower  curve  was  taken  with  a bias  voltage  of 
5.  0 kV.  However,  the  two  curves  have  very  nearly  the  same  slope,  and  should  es- 
sentially lay  on  top  of  one  another  if  they  had  both  been  taken  at  the  same  bias  voltage. 

There  is  a departure  from  this  trend  at  very  early  injection  times  on  the  upper  curve.  l] 
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Figure  7-42. 


Voltage,  lanear  Injector. 


Voltage,  Linear  Injector. 


However,  this  is  probably  at  least  partially  due  to  the  fact  that  the  magnetic  field  is 
changing  rapidly  early  in  time  so  that  the  average  magnetic  field  during  the  decay 
period  is  somewhat  higher  than  the  plotted  points.  The  data  of  Figure  7-44  is  inter- 
esting since  it  shows  an  approximately  linear  scaling  of  containment  time  with  mag- 

2 

netic  field.  The  containment  time  would  be  expected  to  scale  with  B unless  magnetic 
field  errors  or  anomalous  mechanisms  are  dominating  the  loss  processes.  It  is, 
therefore,  likely  that  the  magnetic  field  errors  play  an  important  part  in  the  electron 
loss  processes  for  this  injector  and  vacuum  vessel  geometry. 

Figure  7-45  shows  the  variation  of  the  potential  decay  with  the  width  of  the 
injection  wdndow.  The  decay  rate  is  found  to  improve  with  shorter  injection  ulndows. 

This  is  probably  due  to  the  fact  that  some  neutral  gas  is  liberated  during  injection  by 
electron  bombardment  of  the  injector  and  nearby  vacuum  wall.  Reducing  the  width 
of  the  injection  window  reduces  the  amount  of  gas  liberated  and  hence  the  ion  com- 
ponent in  the  electron  column,  resulting  in  longer  containment  time. 

In  Figure  7-46  the  effect  which  varying  the  pressure  has  on  the  decay  rate  are 
showm.  This  data  was  obtained  by  admitting  small  amounts  of  methane  gas  into  the  chamber 

prior  to  firing  the  machine.  As  can  be  seen,  a systematic  improvement  in  the  con- 

“"G 

tainment  time  is  obtained  as  the  pressure  is  varied  from  1x10  Torr  down  to 
-9 

8x10  Torr.  Figure  7-47  shows  the  effects  of  pressure  variation  on  the  peak  poten- 
tial well  depths  achieved.  The  results  show  a marked  decrease  in  the  potential  well 
depth  at  higher  pressures. 

Some  preliminary  beam  acceleration  experiments  were  performed  with  the 
low  voltage  linear  injector  to  investigate  the  physics  of  the  acceleration  mechanism. 

Figure  7-48  shows  oscilloscope  traces  of  the  beam  acceleration  mode  using  the  low 
voltage  linear  injector.  The  upper  trace  t.hows  the  magnetic  field,  the  second  trace 
shows  the  accelerating  electric  field,  and  the  lower  trace  shows  the  beam  current 
measured  with  a Rogowski  loop  on  a sweep  speed  of  20  p sec/division.  The  injection 
phase  was  started  2.2  msec  after  the  start  of  the  magnetic  field.  The  injection  window 
width  was  —100  psec,  and  the  acceleration  phase  was  begun  at  2300  psec  which  corres- 
ponds to  the  end  of  the  injection  window.  The  beam  current  pulse  width  of  ~30  psec 
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Figure  7-46.  Variation  of  Decay  Rate  with  Pressure  for  the  Low 
Linear  Injector. 


Figure  7-47.  Potential  Well  Depth  as  a Function  of  Pressure  for  the  Low 
Voltage,  Linear  Injector. 
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is  consistent  \\ith  the  time  to  accelerate  the  electrons  from  rest  to  y ~1. 05  to  1. 00. 

At  this  value  of  y the  electron  drift  to  the  wall  due  to  the  toroidal  magnetic  field 
gradient  is  more  rapid  than  the  precession  around  the  minor  axis  which  cancels  the 
drift  effects.  As  a result,  the  electrons  are  lost  to  the  wall.  That  is,  more  charge 
is  required  to  proNdde  the  rotational  transform  needed  for  higher  y.  This  point  is 
further  illustrated  in  Figure  7-49  which  shows  the  theoretical  and  measured  current 
pulse  width  as  a function  of  the  accelerating  electric  field  E ,.  As  can  be  seen,  the 

V 

current  pulse  widths  are  consistent  with  y ~1.05  to  l.OG.  This  agrees  with  theoreti- 
cal predictions  of  the  maximum  y which  can  be  supported  by  ~12  picoulomb  charge. 

The  accelerated  beam  appears  to  be  stable  and  well  behaved  for  the  duration  of  the  20 
to  30  (isec  pulse  width.  However,  the  peak  amplitude  of  the  current  pulse  is  consider- 
ably lower  than  would  be  expected.  For  example,  a total  of  12  ^coulombs  of  charge  ac- 
celerated to  y~  1.05  corresponds  to  a current  of  ~350  amps,  and  only  -50  amps  was 
observed.  This  indicates  that  only  part  of  the  charge  is  being  accelerated.  This  could 
be  due  to  magnetic  field  inhomogenieties  resulting  in  the  trapping  of  the  charge  in  local 
mirrors. 

In  summary,  using  the  low  voltage,  linear  injector,  electron  densities  of 
9 *"3 

~10  cm  were  achieved.  This  corresponds  to  an  injected  charge  of  ~12  pcoulombs 
and  a potential  well  depth  of  ~33  kV.  Charge  containment  times  (time  to  1 'e)  of  sev- 
eral hundred  psec  were  observed,  and  preliminary  beam  acceleration  experiments 
showed  current  pulse  widths  in  agreement  with  theory,  but  beam  currents  which  were 
anomalously  low.  The  injected  charge  w'as  found  to  scale  linearly  with  the  injector 
bias  voltage  for  the  linear  injector.  Based  on  these  results,  modifications  were  made 
to  the  toroidal  vacuum  vessel  to  improve  the  magnetic  field  uniformity,  and  at  the  same 
time  modifications  were  incorporated  to  allow  for  the  use  of  higher  voltage  injectors 
in  order  to  increase  the  injected  charge. 

7.2.4  Vacuum  Vessel  Modifications 

Figure  7-50  shows  a comparison  of  the  STP  vacuum  vessel  configuration  before 
and  after  modification.  For  clarity,  some  of  the  ports  not  involved  in  the  modification 
have  been  left  off  of  Figure  7-50.  The  modification  consisted  of  removing  the  ceramic 
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sector,  the  bellows,  the  filament  sector,  and  the  three  aluminum  "O"  ring  flanges. 

The  ceramic  sector,  the  bellows  and  the  filament  sector  were  replaced  with  smooth 

toroidal  sections.  In  place  of  the  three  aluminum  "O"  ring  flanges,  two  flanges  were 

added  which  provide  electrical  insulation  around  the  toroidal  direction  and  in  addition 

provide  a much  more  uniform  cross  section  than  the  original  flanges.  A sketch  of 

these  flanges  is  showTi  in  Figure  7-51.  This  flange  uses  a viton  "O"  ring  and  .035 

inch  thick  G-10  spacers  to  provide  the  electrical  insulation  provided  by  the  ceramic 

sector  in  the  original  configuration.  The  flanges  are  held  together  by  two  rotatable 

G-10  rings  which  press  against  six  0.15  inch  high,  3/4  inch  wide  lugs.  The  two  rings 

are  held  together  by  3/8  inch  diameter  permali  bolts.  These  flanges  result  in  almost 

no  effective  increase  in  the  amount  of  metal  at  the  flange  location  compared  to  the 

torus  wall  itself  and,  as  a result,  should  lead  to  no  magnetic  field  perturbations. 

The  use  of  viton  "O"  rings  has  not  significantly  degraded  the  operating  vacuum  which 

-9 

is  typically  in  the  low  10  Torr  range.  In  addition  to  these  modifications,  two  4 inch 
OD  ports  were  added  as  shown  in  Figure  7-50b.  These  ports  are  designed  to  give 
ample  clearance  for  high  voltage  injectors  of  the  type  which  will  be  described  in  the 
following  section. 

7.2.5  High  Voltage  Injector  Development 

By  experimenting  with  a number  of  different  geometries,  an  injector  was  de- 
veloped which  injects  more  than  100  pcoulombs  of  charge  into  the  STP  machine.  Some 
of  the  injector  geometries  studied  are  shown  in  Figure  7-52,  and  the  results  of  charge 
injection  measurements  are  given  in  Figure  7-53.  In  addition,  the  effects  of  dimen- 
sional changes  were  investigated  for  some  of  the  injector  geometries  shown  in 
Figure  7-52.  As  can  be  seen  from  Figure  7-53,  the  configuration  of  Figure  7-52(d) 
gave  the  largest  injected  charge. 

On  the  basis  of  these  tests,  a high  voltage  injector  was  built  using  the  geometry’ 
on  Figure  7-52(d).  A photograph  of  this  device  is  shown  in  Figure  7-54,  and  a sketch 
is  shown  in  Figure  7-55.  This  device  uses  a four  conductor  high  voltage  vacuum  feed 
through  to  bring  the  injector  voltages  and  currents  into  the  machine.  These  feed 
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Figure  7-52.  Side  Views  of  V'arious  Injector  Geometries  Used  in 
Charge  Injection  Studies. 
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Figure  7-53.  Summary  of  Results  of  Injector  Geometrj’  Study. 


Figure  7-54.  Photograph  of  the  STP  High  Voltage  Injector, 


THREADED  ROD  FOR 
BELLOWS  ADJUSTMENT 


Figure  7-55,  Sketch  of  High  X'oltage  Iniector. 


throughs  are  mounted  on  a 6"  OD  vacuum  flange  and  the  entire  assembly  is  designed 
to  fit  inside  a 4 inch  OD  port.  This  large  port  size  allows  large  distances  between 
the  high  voltage  injector  components  and  the  torus  walls  which  are  at  ground  potential. 
The  feed  through  flange  is  mounted  on  a bellows  which  is  mounted  on  a vacuum  flange 
at  the  machine.  The  bellows  allows  the  injector  position  in  the  magnetic  field  to  be 
changed  without  disturbing  the  vacuum.  The  injector  itself  consists  of  a tungsten 
dispenser  cathode  and  two  plates.  The  dispenser  cathode  is  a seven  turn,  1 /4  inch 
OD  helical  filament  whose  axis  is  parallel  to  the  toroidal  minor  axis.  The  injector 
voltage  is  applied  between  the  plates  and  the  dispenser  cathode  giving  a configuration 
where  the  electric  field  due  to  the  injector  voltage  is  perpendicular  to  the  toroidal 
magnetic  field. 

7.2.6  High  Voltage  Injector  Results 

Figure  7-56  shows  potential  probe  measurements  of  the  potential  well  depth 
due  to  charge  injected  with  the  high  voltage  injector.  This  data  shows  that  the  poten- 
tial well  depth  and  hence  the  injected  charge  increases  linearly  with  injector  bias  volt- 
age. Potential  wells  of  ~ 300  kV  were  achieved  with  ~35  kV  injector  voltages.  The 
voltage  in  the  center  of  the  machine  could  not  be  measured  directly  above  ~100  kV  be- 
cause the  potential  probes  tended  to  flash  over  between  80  and  100  kV.  However,  the 
voltage  at  the  center  of  the  machine  can  be  extrapolated  from  measurements  made  near 
the  wall.  In  order  to  validate  this  technique  a number  of  radial  profiles  were  taken  at 
various  lower  injector  voltages  to  ensure  that  the  profile  did  not  change  appreciably 
as  the  bias  voltage  changed.  A typical  radial  potential  profile  taken  with  the  high  volt- 
age injector  is  shown  in  Figure  7-57. 

The  charge  injected  with  the  high  voltage  injector  is  plotted  in  Figure  7-58  as 
a function  of  the  injector  bias  voltage.  As  can  be  seen,  in  excess  of  100  p coulombs 
was  obtained  with  ~35  kV  on  the  injector.  In  some  cases  the  charge  deduced  from 
diocotron  frequency  measurements  was  somewhat  lower  than  that  taken  from  potential 
probe  measurements.  However,  the  diocotron  frequency  measurements  were  made 
late  In  time  since  the  diocotron  oscillations  generally  did  not  start  until  several  hundred 


136 


microseconds  after  the  end  of  the  injection  phase.  As  a result,  it  was  necessary  to 
extrapolate  the  density  measured  at  this  time  back  to  the  peak  density  using  the  poten- 
tial probe  decay  rate.  This  limits  the  accuracy  of  the  diocotron  measurements.  Con- 
sidering this,  the  agreement  obtained  between  the  two  measurements  is  reasonable. 

In  addition  to  the  measurements  described  above,  measurements  were  also 
made  of  the  microwave  radiation  fi-om  the  machine.  These  were  made  using  \VR-90 
x-band  waveguide  (~1"  x 0.5")  which  has  a lowest  order  mode  (TEj^^)  pass  band  of 
8.2  to  12.4  GHz.  Figure  7-59  shows  the  microwave  circuit  used.  The  waveguide  ex- 
tended from  a glass  vacuum  view  port  mounted  on  the  torus  to  the  screen  room  some 
50  feet  away.  A variable  attenuator,  ferrite  isolator,  and  diode  detector  were  located 
inside  the  screen  room. 

Figure  7-60  shows  oscilloscope  traces  of  detected  microwave  signals  obtained 
from  the  experiment  for  two  different  injected  charge  levels.  The  upper  photograph 
shows  the  microwave  signal  for  a case  where  the  injected  charge  was  ~42  pcoulombs 
while  the  lower  photograph  corresponds  to  ~81  p coulombs  of  injected  charge.  As  can 
be  seen  from  the  photographs,  both  the  amplitude  and  shape  of  the  microwave  signal 
changes  as  the  electron  density  changes.  In  both  cases,  the  microwave  emission  is 
strong  at  the  beginning  of  injection.  However,  at  the  lower  charge  level  the  emission 
dies  away  rather  rapidly  after  the  early  peak,  while  at  higher  injected  charge  levels 
the  emission  has  a longer  duration.  Since  the  emission  is  strongest  in  both  cases 
early  in  time  when  the  electron  density  has  not  reached  its  peak,  we  assume  the  early 
peak  is  due  to  a collective  radiation  mechanism  driven  by  the  injection  process.  The 
emission  occuring  later  in  time  is  probably  due  to  ordinary  incoherent  cyclotron  radi- 
ation from  the  electrons  in  the  toroidal  magnetic  field.  The  electron  density  at  the 
higher  charge  levels  is  probably  sufficient  to  give  measurable  incoherent  cyclotron 
emission  while  at  the  lower  densities  the  radiated  power  through  the  small  vacuum 
port  is  probably  too  small  to  be  detected. 

The  frequency  of  the  radiated  microwave  emission  should  depend  on  the  strength 
of  the  magnetic  field  if  the  emission  occurs  at  the  cyclotron  frequency.  The  electron 


Figure  7-()0,  Oscilloscope  I'races  of  Radiated  Microwave  Power  From  the  S'!'!’ 

Machine,  i'he  upper  trace  in  each  photograph  is  the  microwave 
signal,  and  the  lower  trace  in  each  case  is  a potential  probe 
signal.  All  traces  are  30  ixsec  div.  The  upix^r  microwave  signal 
sensitivity  is  2 mv/div  and  the  lower  is  3 mv/div.  The  upper 
photograph  corresponds  to  ~12  jufoulombs  of  injected  change  while 
the  lower  corresponds  to  ~81  jit'oulombs. 


cyclotron  frequency  is  given  by 


This  corresponds  to  ~2. 8 GHz  per  kilogauss  of  magnetic  field.  Thus,  by  changing 
the  toroidal  magnetic  field  in  the  experiment  it  is  possible  to  vary  the  cyclotron  fre- 
quency through  the  x-band  range.  Figure  7-61  shows  a plot  of  the  detected  microwave 
signal  as  a function  of  the  cyclotron  frequency  calculated  for  an  electron  at  the  minor 
axis  of  the  machine  as  the  magnetic  field  is  varied.  The  magnetic  field  was  varied 
both  by  changing  the  peak  magnetic  field  and  by  changing  the  injection  time,  giving 
nearly  identical  results.  The  fact  that  the  detected  microwave  power  depends  on  the 
magnetic  field  indicates  that  the  emission  is  probably  cyclotron  radiation. 

7.2.7  Charge  Acceleration  Studies 

The  last  few'  weeks  of  this  year's  work  were  devoted  to  charge  acceleration 

experiments.  These  efforts  indicated  some  unanticipated  problems  in  the  charge 

acceleration  phase.  Although  ample  charge  was  available  to  provide  the  E and  B 
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particle  rotational  transform  to  cancel  the  toroidal  drift,  and  the  vertical  magnetic 
field  was  available  for  beam  |X)sitioning,  no  large  un-neutralized  beam  curreiils  were 
established.  Currents  initially  thought  to  be  electron  beam  currents  were  later  found 
to  be  flowing  in  the  wall  of  the  vacuum  vessel,  by-passing  one  of  the  insulating  vacuum 
flanges  which  was  slightly  misaligned.  In  some  cases,  large  currents  were  observed 
late  in  time  due  to  ionization  of  wall  material  and  background  gas  to  form  a plasma 
allowing  a tokamak  discharge.  Figure  7-62  shows  oscilloscope  traces  of  a tokamak 
discharge.  As  can  be  seen,  a discharge  current  of  ~10  kA  peak  occurred  late  in 
time  (~2msec  after  injection).  The  x-rays  generated  near  the  peak  of  the  current 
pulse  were  found  to  penetrate  up  to  ~0. 25  inches  of  lead,  indicating  the  presence 
of  electrons  with  up  to  ~200  keV  energy  in  the  discharge.  This  implies  that  these 
electrons  were  contained  for  several  thousand  passes  since  the  transformer  electric 
field  was  ~0. 2 volts/cm  which  corresponds  to  an  energy  increase  of  ~64  eV  per 
transit  for  the  electrons.  The  long  delay  before  the  onset  of  the  tokamak  discharge 
is  probably  the  result  of  a slow  buildup  of  plasma  density  by  ionization  of  background 
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Figure  7-62.  Oscilloscope  Traces  of  a Late  Time  Tokamak  Discharge  in 

the  STP  Machine.  Upper  trace:  Toroidal  current,  14  kA/div. 
Lower  trace:  Scintillator  - photomultiplier  x-ray  signal. 
Horizontal  sensitivity:  500  usec/div.  The  oscilloscope  was 
triggered  at  injection  time.  The  transformer  electric  field 
was  ~0. 2 volts/cm. 


gas  and  wall  material.  Normally,  one  of  the  insulating  torus  flanges  (see  Figure  7-51) 
was  short  circuited  by  ground  connections  between  the  two  halves  of  the  torus,  leaving 
only  one  insulating  gap  in  the  vacuum  vessel.  Under  these  circumstances  the  tokamak 
discharges  were  usually  observed  late  in  time.  When  the  ground  was  removed  from 
half  of  the  machine,  thus  removing  the  shorting  path  around  one  of  the  flanges,  no  toka- 
mak discharge  occurred.  This  indicates  that  one  of  the  flanges  drew  some  leakage 
current  generating  enough  plasma  for  the  tokamak  discharge.  This  leakage  current 
was  observed  to  begin  early  in  time  and  persist  until  the  onset  of  the  tokamak  mode. 

The  amplitude  of  this  leakage  current  was  a few  hundred  amps. 

From  the  results  of  these  experiments  we  conclude  that  the  machine  operates 

as  expected  in  the  presence  of  a plasma,  generating  a tokamak  discharge,  but  does 

not  generate  unneutralized  electron  beam  currents  even  when  a large  amount  of  charge 

is  present.  (The  charge  decay  time  is  not  affected  by  the  presence  of  the  toroidal 

electric  field. ) This  suggests  that  the  difficulties  are  connected  with  the  orbits  of  the 

electrons  in  the  unneutralized  column.  These  orbits  involve  principally  motion  around 

the  toroidal  minor  axis  (due  to  E x B drift)  so  that  the  particles  can  be  expected  to  have 

o 

much  greater  energy  perjjendicular  to  the  magnetic  field  than  parallel  to  the  field.  Such 
particles  would  be  susceptible  to  being  trapped  in  small  ripples  in  the  magnetic  field 
caused  by  the  discrete  nature  of  the  magnetic  field  coils.  Most  of  the  electrons  trapped 
in  such  wells  would  not  feel  the  effects  of  the  toroidal  electric  field  since  the  electric 
field  would  result  only  in  a small  current  even  though  a large  amount  of  charge  is  pre- 
sent for  several  hundred  microseconds  after  injection.  This  is  precisely  the  behavior 
that  is  observed.  The  currents  observed  are  ~50  amps. 

This  situation  can  be  corrected  by  injecting  electrons  with  significant  energy 
parallel  to  the  magnetic  field  that  they  will  not  be  trapped  in  magnetic  wells.  Such 
electrons  would  then  be  further  accelerated  by  the  transformer  each  time  they  passed 
under  the  insulating  vacuum  flanges.  This  injection  process  could  be  best  accomplished 
with  a relativistic  Injector.  Figure  7-63  shows  a drawing  of  a relativistic  electron  in- 
jector. Basically  the  device  would  consist  of  a diode  driven  by  a marx  generator.  The 
diode  is  connected  to  the  marx  by  a vacuum  coaxial  line.  The  insulator  and  marx 
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Insulator  and  Vacuum  Coax  Layout  for  the  STP  Machine. 


generator  would  be  located  in  an  oil  tank,  and  a rough  vacuum  would  be  used  as  a 
buffer  region  between  the  oil  seal  and  the  high  vacuum  region.  The  marx  generator, 
oil  tank,  and  insulator  for  this  injector  are  part  of  an  existing  small  generator  which 
can  be  readfly  adapted  to  relativistic  injection.  The  marx  would  be  equipped  with  a 
crowbar  so  that  the  voltage  could  be  removed  from  the  diode  before  gap  closure 
occurred. 

Computer  calculations  indicate  that  this  scheme  could  lead  to  a large  increase 
in  the  charge  level  inside  the  machine.  A similar  injection  scheme  has  been  used  by 
researchers  in  Japan^  to  produce  an  unneutralized  electron  beam  in  a tokamak-like 
device.  With  the  presence  of  a previously  injected  thermionic  core  of  charge  the  in- 
jection and  trapping  efficiency  in  the  STP  machine  should  be  good.  Thus,  this  approach 
should  lead  to  both  higher  charge  levels  and  significant  beam  currents. 


^ A.  Mohri,  M.  Masuzake,  T.  Tsuzuke,  and  K.  Ikuta.  "Formation  of  a Non-Neutral 
Relativistic  Electron  Beam  Ring  in  a Toroidal  Magnetic  Field,  Phys.  Rev.  Letters, 
34,  574  (1975). 
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